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HAUPTAUFSÄTZE 


Über die Potentialströmungen zäher Flüssigkeiten. 


Abhandlungen zur Hydrodynamik IX. 
Von Georg Hamel in Berlin. 


Es werden die Potentialströmungen zäher Flüssigkeiten. die es bekanntlich 
gibt, die aber des Haftens am Rande wegen gewöhnlich beiseite gelassen 
werden, genauer auf den Energieaustausch hin untersucht. Für die ebene 
Bewegung findet sich die einfache Formel (9) der inneren Leistung, die 
ebenso einfache Formel (10) der äußeren Leistung. Natürlich sind beide 
entgegengesetzt gleich, was aus dem komplex geschriebenen Greenschen 
Satz sofort folgt. Die Randbedingung des Haftens muß durch Nachgeben 
des eingebetteten Körpers erreicht werden. Deshalb erscheint die im zweiten 
Teil behandelte Frage von Interesse, ob es Zylinder gibt mit konstanter 
Winkelgeschwindigkeit am Rande. Dieses Problem führt auf eine quadra- 
tische Randwertaufgabe (4, (5), von der aber nur die triviale Lösung des 
Kreiszylinders angegeben werden kann. Nachbarlösungen dazu gibt es in 
einem bestimmten Sinne nicht. 


Bekanntlich sind Potentialbewegungen volumbeständiger Flüssigkeiten mögliche Be- 
wegungen zäher Flüssigkeiten; sie werden aber meist außer acht gelassen, da erstens zähe 
Flüssigkeiten an den Wänden haften müssen, daher nach bekannten Sätzen der Potential- 
theorie die Flüssigkeit in Ruhe bleibt, falls die Wände sich sicht bewegen und da zweitens 
die im Falle der Potentialbewegung auch für zähe Flüssigkeiten geltende Bernoullische 
Gleichung zu zeigen scheint, daß ein Energieaustausch nicht anders stattfindet als bei idealen 
Flüssigkeiten. 


Beides ist ein Vorurteil: erstens gibt es auch bewegte Wände, die nicht starr zu sein 
brauchen, sich also umgekehrt der Bewegung der Flüssigkeit anpassen können, und zweitens 
gibt die Bernoullische Gleichung nieht den wahren Energieaustausch, tut dies vielmehr 
nur für eine Ersatzdynamik, in der die wahren Kräfte durch ihre Resultierende ersetzt sind. 
Und das ist energetisch nicht dasselbe, wie elementare Beispiele leicht zeigen. 


Im folgenden soll zuerst der Energieaustausch, namentlich im Falle der ebenen Potential- 
bewegung zäher Flüssigkeiten genauer untersucht und an Beispielen erörtert werden. Zweitens 
stellen wir die Frage, ob es zylindrische Körper gibt, um die eine Potentialbewegung so statt- 
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findet, daß an der Oberfläche die Winkelgeschwindigkeit konstant ist. Diese Frage führt auf 
eine neuartige Randwertaufgabe der Funktionentheorie. Es zeigt sich, daß unter gewissen 
Voraussetzungen der Kreiszylinder die einzige Lösung ist, wenigstens in dem Sinne, daß es 
zu ihm keine Nachbarlösungen gibt. (Die genaue Formulierung des Satzes am Schlusse der 
Arbeit.) 


Teil I. Der Energieausgleich. 


1. Bekannte Tatsachen. Die Bewegungsgleichungen zäher Flüssigkeiten lassen sich be- 
kanntlich so schreiben: 


Od 
dt 
wo b den Vektor der Geschwindigkeit, p den Druck, o die konstante Dichte, » die (konstante) 


kinematische Zähigkeit, % die Höhe, g die Erdbeschleunigung bedeutet. Dazu kommt die so- 
genannte Kontinuitätsgleichung 


Man erkennt sofort die Möglichkeit der Potentialbewegung mit 
v—grady und also rotv=0, 
wodurch sich (1) zur Bernoullischen N 


dt 
integriert. Denn es ist ja Iv= Igrady=grad Jy=0, was wegen (2) zutrifft. 


Das Glied » Av rührt von den inneren Spannungen der Zähigkeit her, von dem Tensor 
%, = 2, de... Z,, wobei 


n ist der Koeffizient der inneren Reibung, „=or. v,, t,, v, sind die Koordinaten des 


Vektors d. Nun sind bekanntlich die Spannungskoordinaten X,, Y„, Z„ für ein Flächen- 
element d F mit der äußeren Normalen n 


(n, x) + X,cos (n, y)-+ X. cos (n,z) usw. 


und somit ist die von außen auf einen abgeschlossenen Teil der Flüssigkeit übertragene 
Leistung der Zähigkeitsspannungen 


Die innere Leistung ist aber bekanntlich das Volumintegral 


und das ist stets . ee nur Null bei starrer Bewegung. In allen anderen Fällen wird 
also von der inneren Reibung im Innern der Flüssigkeit Energie verbraucht. 
Wenn nun ein Potential g da ist, nimmt (6) die Form an 


2. Leistungsformeln im Falle der ebenen Potentialbewegung. Im r'alle der ebenen Be- 


wegung ist 


1) Siehe z. B. Lamb-Friedel: Lehrbuch der Hydrodynamik. Leipzig und Berlin (1907), S. 668. 
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und wenn wir das komplexe Potential 


y+tiy=fle+tiy 
einführen, ist 


| 


Daher ist 
und wir bekommen für die innere Leistung die einfache Formel » 


Alle diese Integrale sind hier und hinfort je Einheit der dritten Dimension verstanden. 
Auch für die äußere Leistung (5) bekommen wir einen überaus einfachen Ausdruck. 
Aus (5) wird zunächst mit ds als Bogenelement 


+ + X,ecos(n, y)] 
+[Yx cos + Y„eos (n,y)]e,y ds + (Xyr. + Y,r)ds 


dyt 


- ut 


-f'fidz, 
was wir auch formulieren können: Die äußere Leistung der Zähigkeit ist 
Man kann diese Formel noch umformen. Es sei: 
=v,-ivy=verid, 
df’=e-id(de -ivdd), 
faf- faf=r(de+iedd) e(de -irdd)=2indd 
und also 
L=- v dd 
Das scheinbare Herausfallen der Zähigkeitswirkung in der Bernoullischen Gleichung findet 
hier ihr Gegenstück in dem Satz, daß bei Potentialbewegung 
ist. Denn aus (5) folgt 
La= Yuvy+ Zuv)dF= 0x + eos (n,a) +: 


\ 12 (Ar + + v,) + dy (X, + Y, + Zy v.) + oz (X, Ux + % "y + Z, 


/ 


?2) Lamb gibt S.669 (1 bis 6) für den Raum die Formel — 3 "Bar. 
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dx 
weiteren Glieder geben — L;. Die Sache ist nach u allgemeinen Prinzipien der Mechanik 
überdies klar. Ferner ist der Satz nichts anderes, als der in Komplexe umgeschriebene Green- 


sche Satz: as. 


Die drei ersten Glieder verschwinden, da -—=ndv,—0 ist usw., und die 


3. Beispiele. Als Beispiele nehmen wir 
1. Die Drehbewegung um einen Kreiszylinder. Im Unendlichen soll Rube 
herrschen. Dann ist 


(# +iy), € reell konstant, 


die einzige Lösung; sie ist auch die einzige überhaupt, wenn wir alle Lösungen von Rotations- 
symmetrie hefanziehen, nämlich 


v=Ar’nr+Br+Chr+D, 


von denen die beiden ersten Glieder wegen der Bedingung im Unendlichen fortfallen müssen. 
Es ist also 


@_ (dr 
a 


ci 
Die Umfangsgeschwindigkeit ist also 


Dagegen ist über einen Kreis integriert 


Dies verschwindet für r>»x, von außen wird also keine Leistung übertragen, 
dies geschieht allein von dem sich drehenden Zylinder, der die ganze ver- 
brauchte Leistung 477 aufbringen muß. 


2. Als zweites Beispiel nehmen wir den umströmten Kreiszylinder; die 
Anströmgeschwindigkeit im Unendlichen sei jetzt U, sie erfolge in Richtung der x»- Achse. 
Dann ist bekanntlich die gie gegeben durch 

2 


f= (e+iy®' 


Also die innere Leistung 


a 


Diese Leistung muß auch hier lediglich durch den Zylinder aufgebracht werden. Denn die 
äußere Leistung an einer Zylinderfläche ist 


4 4 
= 422 07; 


Dieses gibt für r>» wieder null: die äußere Strömung gibt die Leistung nicht 
her, der Zylinder muß die ganze Leistung 8xnÜ? liefern. 
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Dieser Zylinder kann natürlich nicht starr sein; er muß die Bewegung 
der Flüssigkeit mitmachen. Diese ist an seiner Oberfläche durch 


also 
= U(l—cos2d9), vy,=— Usin29 


gegeben. Man müßte also, um die Bewegung zu realisieren, der Haut des Zylinders eine Be- 

wegung erteilen, wie sie Bild 1 andeutet und die entfernt an die Bewegung eines Schwimmers 
erinnert. Kehren wir die Bewegung um und lassen den Zylinder 
sich durch die Flüssigkeit, die jetzt im Unendlichen ruhen möge, 
mit der Geschwindigkeit U bewegen, dabei die angedeutete 
„Schwimmbewegung* ausführen, so müßte er, solange die lami- 
nare Bewegung stabil ist, einen Widerstand von der Größe 877 U 
finden, da die Leistung 87 7 U? die obige bleibt. Größenordnungs- 
gemäß stimmt sie mit einer Formel von Öseen und Lamb über- 
ein, die für den starren Zylinder den Widerstand 


SanU 
1,7811 Re 


Bild ı. 1-2log 


angibt (Re die Reynoldssche Zahl). Beide Formeln würden identisch sein für R, = i rn 

also etwas mehr als Re=4, während für die Gültigkeit der Formel Re= 1 angegeben wird. 
3. Als drittes Beispiel nehmen wir die Strömung um zwei gleiche 

parallele Kreiszylinder; im Unendlichen herrsche Ruhe. 2c sei die Entfernung 

(ler Mittelpunkte voneinander; R sei ihr Radius. Sie 

entstehen durch die konforme Abbildung 

‚_2—a 


aus zwei konzentrischen Kreiser 


1 
und-Iw 


der w-Ebene; man berechnet leicht e=Aa; R=ay#? -1; 
Wir können annehmen: 


entsteht aus o? durch Vertauschen von mit — 4. 


Strömungspotential ist gegeben durch 


a 
a’ 


C C 2 — 


daher 


Daraus ergibt sich die Geschwindigkeitsverteilung an der Oberfläche der Zylinder. 


Im Punkte 0 ist 2=0, also -ir,„= Sri also v.=0, 1, = — 
Wir setzen diesen Wert gleich U, so daß 
C=-aaU=-— U=-aye-RU. 


Die verbrauchte Leistung kann man nun entweder aus 


oder aus 
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berechnen. Man erhält 


\4ra m? mw? 


27 


zn 


we o* odod# (m — ei”) 


= 


L. berechnet sich so: 
ff’da) 


über beide Zylinder erstreckt; oder 


zilz 
La =2n 4a 3 (2? - a?) (z? a?) 
über den rechten Zylinder erstreckt; die Verdoppelung ist durch den ersten Faktor 2 vor- 
genommen worden. 
Nun heißt aber die Kreisgleichung 22 —/a(z2+2)+.«a”—=0. Eliminiert man damit z, 
so bekommt man 


(A? — (2? — a?) (2? - u) 


Dieses Integral berechnet man nach dem Residuensatz zu 


Dieses ist die Leistung, welche die strömende Flüssigkeit an die Zylinder abgibt; daher das 
- Zeichen. 


4. Als viertes Beispiel nehmen wir die Strömung um ein Joukowski- 
Kutta-Profil. Es liege in der 2-Ebene und werde auf den Kreis 'Ü{|=R in einer {-Ebene 
durch die bekannten Transformationen 

(mw w,); w=2 + 2’; w,— 
konform abgebildet. 0 in der w-Ebene ist der Punkt w= 0; w=R, d.h. 
z=R, Re’? entspricht der hinteren Spitze des Tragflügels; w—= — R, 
d.h. R und 


dem durch Abrundung verdeckten vorderen Ende des Kreisbogens, dessen 
Bild ein den gezeichneten Kreis in »= R berührender Kreis ist. 
Sind nun die Anströmkoordinaten U und V, so wäre für den Kreis Bild 3 und 4. 
der Ö-Ebene das Potential bekanntlich 


Da aber für große { resp. z der Grenzwert von 3 gleich > ist, muß man, wenn U- iV die 


Anströmgeschwindigkeit in der z-Ebene bleiben soll, mit SR multiplizieren. Somit erhält man 
für die z-Ebene die komplexe Geschwindigkeit 
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R 


= -iN-(U+iV) 
Bestimmt man nun in bekannter Weise € so, daß die Singularität am hinteren Ende w — R 
fortfällt, d. h. setzt man 
C=2ai 
so bekommt man nach leichter Rechnung 


- - 


Bemerkenswert ist, daß dies eine rationale Funktion von £ ist, deren Nenner Z?(Ö  [,) ist, 


deren Zähler ebenfalls dritten Grades ist, so daß = die Partialbruchzerlegung 


hat. A brauchen wir nicht, es berechnen sich 


wi (U+iV) Re, 
D= 


Nun kann man die Formel für die äußere Leistung leichtkomplex trans- 
formieren. Es ist 


Wir können also das Integral in der {-Ebene über den Kreis |{| = R erstrecken, wo dann 


2 


fa DR ist. f’ ist unverändert das obige 2 also 


dz 


oder auf dem Kreis !Ü\= R 
R 


2 
Der Pol von f’, nämlich £ -? liegt außerhalb. 


Mithin = der Residuensatz 
BB „co BD BB CD, DDR: 
$r‘ f' 2 2 zur 2 


R R® R® (R?— 5) 
Also ist das — 4 77-fache der vorstehenden eekigen Klammer. 


2 


R 
Wir wollen nun den echten Bruch =. setzen. Dann ist der Nenner 


27 Re") = 4 R’e(cosp — e). 
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g g 


Die positive Größe cos (acos$— R) hat ersichtlich eine besondere Be- 


deutung für den Tragflügel; sie ist ein Maß für die Abrundung. 
Eine längere, aber elementare Rechnung ergibt nun die Endformel: 


e? cos — V?]- (cos ß— | > 
Da die Abrundung in der Regel gering sein wird, ist cosß — e klein, also das letzte 
Glied das größte. Behält man es allein bei und setzt man im Zähler überall e= cos p, so 
bekommt man als Näherungsformel für wenig abgerundete Tragflügel 
(U cos + V sin p)? + cos?  V? 


Führt man noch den Anstellwinkel « durch UV=»cosa, V=rsina ein, so bekommt man 
_ — A) + cos? Psin?a 

La=-- ran v? cos?’ (cos 

und also einen Widerstand 

_ — P) + cos? a 

W=anv cos’ (cos 


Es ist offensichtlich nur das Glied berücksichtigt, das von dem scharfen Umbiegen um die 
Vorderkante herrührt. Ein Widerstand, der von Wirbelablösung herrührt, steckt natürlich 
nicht in der Formel; daher ist ein Vergleich mit der Erfahrung kaum möglich, zumal diese 
ÖOberflächenbewegung wohl noch schwerer zu realisieren ist, als die beim Kreiszylinder in 
Beispiel (2). Vielleicht daß das Ergebnis in Analogie mit diesem zweiten Beispiel .eine ge- 
wisse Abschätzung des Widerstandes bei laminarer Bewegung bedeutet. 


Teil I. Zylinder mit konstanter Winkelgeschwindigkeit an der Obeıfläche. 


1. Der Ansatz. Es hat vielleicht Interesse, nach solchen Profilen zu fragen, die mit 
konstanter Winkelgeschwindigkeit durchlaufen werden. Man könnte solche Bewegungen tech- 
nisch leichter ausführen. 

Diese Aufgabe läßt sich auf ein Randwertproblem zurückführen, das von der vierten 
Ordnung ist. Eine Lösung ist bekannt: der Kreiszylinder, unser erstes Beispiel. Es scheint 
aber, daß es weitere Lösungen nicht gibt, jedenfalls gibt es in gewissem Sinne keine stetig 
benachbarten Lösungen zum Kreiszylinder. Vielleicht gibt es Lösungen, die aus mehreren 
Zylindern bestehen, also ein mehrfach zusammenhängendes Strömungsgebiet ergeben. Darüber 
aber kann ich nichts sagen. 


Wir wollen das Problem für einen Zylinder ansetzen. Es soll also sein 
e'’ mit d=w t, konstant 
Mithin 
. .. dr id 
datrir e’o, 
folglich 


Bilden wir nun das Äußere des gesuchten Zylinders auf das Äußere des Kreises |w|=o einer 


w-Ebene durch 
z—= F(w) 


konform ab, so entsteht dort die durch 
C 
f log w 
gegebene Strömung. Es ist also 


dz 2rimwF'(m)' 


Auf dem gesuchten Rand, d. h. für w=oe'“ soll das nach (2) gleich 
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dr dz 
sein. Also ist 
dg 1dw 
jer 
oder 2rivomdz 
dz ’dw 


oder dw dd 2riomw' 


Dies gilt für dw=-ioe Also folgt 
C 


Nun ist aber nach (1) 


Daraus 


.d9 F(mw)dw F’(w)d mw 


F(n)da Fim) da 


iw+ Fin) in—=2iR F(m) m. 


Setzen wir dies in (3) ein, so bekommen wir die Formulierung der Randwertaufgabe. 
Es soll für |w!>o eine analytische Funktion 


z=F(m)—= N a,w”. 


gesucht werden, die auf |w|=o die Randbedingung 
dz E dz C 


dw 2 2n0’o 


erfüllt. Man kann dieser Bedingung leicht die Form geben 
„dadr 


2 


wo dr das überstrichene Flächenelement r?d# in der z-Ebene, ds das Bogenelement ist. 


2. Durchführung. Wir sehen zunächst, was wir mit dem Ansatz (4) erreichen können. 
Wir formen die linke Seite der Randbedingung (5) um. Aus 


=— 
dw 


—1 


wofür wir auch Ib, w” b,w” schreiben mit b,—= (vr 1)a,_,. Das aber ist wegen 


gleich 


u 
\ b, w m =, b, bu =) w" sv—-ın b, 
v0 u | 


(b mit negativem Index sind null zu setzen) 
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re Fin), 
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=, 20 +, 20 


0 
Da der Ausdruck reell ist, können wir ein R davorschreiben und dann unter dem R Zeichen 
en 


in der ersten Summe zum konjugierten Wert übergehen. Setzen wir dann noch w” — an 


so 


erhalten wir 


IX \ 

1 


L 
b, by + 22 w- b, tn 
1 


Mithin hat die Randbedingung die Form angenommen: 


In der geschweiften Klammer steht eine Funktion, die im Außenraum |w| >. regulär 
ist, denn 2 darf dort nirgendwo null sein. Also muß der reelle Teil dieser Funktion dort 


überall null sein, die Funktion selbst also konstant imaginär: 


1 


Daraus folgt durch Integration 


\ 5) 


log 2 log At | Di) log 


\ m n 


v+rn b, 


(A Integrationskonstante) und daraus 


Damit nun z die Form ... (4) hat, muß 
D=0 
und 

sein. Hieraus folgt |b,’= x0o”” und somit die Konvergenz der Reihe ... (4) für |w! >o, falls 
sich nur die b, resp. a, bestimmen lassen. 

Solche Bestimmungsgleichungen bekommen wir nun aus dem Vergleich der beiden Ent- 
wiecklungen für z 


n 

v 


—1 


A=b,=a_, 


> 


4 
v 1 


Awe 
1 v 


Sofort folgt: 


weiter 


Die weiteren Gleichungen sind komplizierter: es sind © viele nichtlineare Gleichun- 
gen mit unendlich vielen Unbekannten 5, und es sind keine Rekursionsformeln. Es 
scheint nicht aussichtsvoll, sie lösen zu wollen. 

3. Nachbarlösungen zum Kreiszylinder. Obwohl man die Näherungsidee auch direkt an- 
wenden kann, schlagen wir einen anderen Weg ein. Es muß 


z=(,Ww 


de?’ 
dw 
Ti - 
wo 
(9) 
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eine Lösung sein. Durch direktes Einsetzen in ... (5) bekommen wir 


270’ 


Wir fragen nun nach Nachbarlösungen, d. h. wir setzen 


> 
(m) 
mit einem hinreichend kleinen e und 9, = w. 
Setzen wir das in die Randbedingung (5) ein, die wir jetzt besser 


dz dz 
dw de ) 
schreiben, so bekommen wir 


Das = (,C, hat sich beiderseits fortgehoben. 
Eine einfache Rechnung liefern nun die Rekursionsdifferentialgleichungen zunächst für 
den Rand 


v+- u +o+r.n 
Rechts laufen die Indizes nur von 1 ab. Daß 9, w ist, ist schon benutzt. 

Nach der Methode der vorhergehenden Betrachtung kann man nun auch die beiden 
ersten Glieder und damit die ganze rechte Seite in die Form R@ (m) bringen, wo G (m) nach 
fallenden Potenzen von m geordnet ist. 

Man kann dann’ weiter beiderseits nach der Schlußweise der vorigen Nummer das WR 
fortlassen und erhält dann wegen D=0 (Formel 6) Differentialgleichungen der Form 


v 


\ 
Y 
+ w w 
0 


mit bekannten ('„,,. Diese läßt sich integrieren zu 


logw 
In = C,„, ow + C,, ‚+ Ca, 2 + 
3 
B,„ ist eine neue 
Da das Glied mit log w unzulässig ist, müssen die Integrationskonstanten so eingerichtet 
werden, daß alle C„,,=0 sind Die ersten Glieder berechnen sich leicht. Unter Beachtung 
von D=0 bekommt man 


B, 


> w 


B, B, B, Bi log 


Da ein Glied mit log w unzulässig ist, muß 


sein. Damit aber fiele das ganze erste Korrekturglied & g, (m) fort, und wenn man nun das e 
als den Kleinheitsfaktor des ersten wirklich vorhandenen Gliedes definiert, entsteht ein 
Widerspruch. 

Es gibt also keine Nachbarlösungen der Form 


n 


0 In (m) 
0 


mit 9, (®)=w, g,(w)=20. Es bleibt dahingestellt, ob es Entwicklungen nach gebrochenen 
Potenzen geben kann. 260 
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Grenzlinien adiabatischer Potentialströmungen. 
Von W. Tollmien in Dresden. 


Durch Beispiele stationärer ebener adiabatischer Potentialströmungen war 
das Auftreten von „@renzlinien“ bekannt, die eine Schar Machscher Wellen 
einhüllen. Unter Zuhilfenahme der Molenbroek-Tschapliginschen 
Transformation werden nunmehr allgemein die Singularitäten von Gasströ- 
mungen an solchen Grenzlinien erörtert. Es wird bewiesen, daß eine Strömung 
in keiner Weise über die Grenzlinie hinaus fortgesetzt werden kann. Daher 
muß die berechnete adiabatische Potentialströmung schon stromaufmwärts 
vor der Grenzlinie zusammenbrechen. 


Bei der Betrachtung stationärer ebener Strömungen eines kompressiblen Gases 
kann man oft Wirbelfreiheit voraussetzen, falls sich das Gas durchweg adiabatisch verhält. 
Freilich ändert die Differentialgleichung für das Geschwindigkeitspotential in Abhängigkeit 
von den cartesischen Ortskoordinaten ihre Art, wenn die Strömungsgeschwindigkeiten in 
manchen Gebieten die Schallgeschwindigkeit überschreiten. Doch bedingt dieser Übergang 
von einer elliptischen zu einer hyperbolischen Differentialgleichung zunächst nur ein ver- 
schiedenes Verhalten der Strömung bei der Fortsetzung kleiner Störungen, wie sie etwa von 
begrenzenden Wänden ausgehen. Im großen und ganzen beeinflußt aber der Durchgang durch 
die Schallgeschwindigkeit die Strömung noch wenig, wie man etwa aus Druckmessungen 
erkennen kann. Dagegen legen viele Beobachtungen den Verdacht nahe, daß eine starke 
Anderung des Verhaltens der Strömung eintreten kann, sobald die Schall- 
geschwindigkeit von der Strömungsgeschwindigkeit erst einmal um ein gewisses Maß über- 
schritten wird. Die in solchen Fällen mehrfach beobachtete schroffe Änderung der Druck- 
verteilung z. B. an Tragflügelprofilen möchte zu dem Schluß verführen, daß die Strömung 
überhaupt nicht mehr wirbelfrei und adiabatisch bleibt. Darauf hat besonders G. I. Taylor') 
hingewiesen. Auf welches Vorkommnis ist nun der Umsturz der Potentialbewegung zurückzu- 
führen, nachdem der Durchgang durch die Schallgeschwindigkeit an sich aus theoretischen und 
experimentellen Gründen als Ursache ausscheidet? Das ist die Frage. 


Nun war ich bei der Konstruktion adiabatischer Potentialströmungen mit einem Über- 
gang aus dem Unterschall- in das Überschallgebiet °) im Jahre 1937 zum ersten Male auf 
eine Erscheinung gestoßen, die das Weiterbestehen einer Potentialströmung zum mindesten 
sehr fraglich machte. In diesen ebenen Strömungen durch gekrümmte Kanäle traten sehr 
häufig nach Überschreiten der Schallgeschwindigkeit Linien auf, an welchen die Strömungen 
verschiedene singuläre Eigenschaften hatten. Diese Linien hatte ich bis zur weiteren Klärung 
des Sachverhaltes vorläufig Grenzlinien genannt, da mir durch sie den adiabatischen 
Potentialströmungen eine natürliche Grenze gesetzt zu sein schien. Das Druckgefälle an den 
Grenzlinien wird unendlich, die Stromlinien haben dort Spitzen. Für besonders wichtig hielt 
ich die Eigenschaft, daß die Grenzlinie nicht nur eine Machsche Welle ®), sondern auch eine 
Enveloppe Machscher Wellen ist. 


Es ist aber keineswegs so, daß nach Überschreiten der Schallgeschwindigkeit unbedingt 
eine Grenzlinie auftritt. Dies sieht man am Beispiel des Potentialwirbels mit konzen- 
trischen, kreisförmigen Stromlinien und derselben Geschwindigkeitsverteilung, wie sie in den. 
Anfängen der Hydrodynamik ohne Berücksichtigung der Kompressibilität gelehrt wird. Um 
ein kreisringförmiges Gebiet mit Überschallgeschwindigkeit ist also das Unterschallgebiet 
gelagert. Auch in meiner angeführten Arbeit findet sich eine Strömung mit Übergang von 
Unterschall- zu Überschallgeschwindigkeiten, ohne daß eine Grenzlinie auftritt (a.a. ©. S. 129). 


F. Ringleb*) hat nun im Jahre 1940 neue Strömungen mit Grenzlinien gefunden. 
Er nannte diese Linien Stoßlinien, was leicht zu Verwechslungen mit den ganz andersartigen 
Verdichtungsstößen führen kann. Zur Berechnung seiner Strömungen benutzte F. Ringleb 


1) G. 1. Taylor: Recent work on the flow of compressible fluids. J. London Mathem. Soc. Bd. 5 (1930), S. 224 
bis 240. Deutsche Übersetzung in Z. angew. Math. Mech. Bd. 10 (1930), S. 334 bis 345. Vgl. auch den Vortrag von 
G. I. Taylor: Weil established problems in high speed flow. Reale Accad. d’Italia, Fondazione Allesandro Volta, Atti 
dei Convegni 5, 30. Sett. — 6. Ott. 1935. Tema: Le alte veloeitä in aviazione S. 198 bis 214, Roma 1936. 


2) W. Tollmien: Zum Übergang von Unterschall- in UÜberschallströmungen. Z. angew. Math. Mech. Bd. 17 
(1937), S. 117 bis 136. 


3) Daß die Grenzlinie eine Machsche Welle ist, hatte schon W.G. Bickley an einem von Taylor herrührenden 
Beispiel einer aufhörenden adiabatischen Potentialströmung bemerkt, vgl Z. angew. Math. Mech. Bd. 10 (1930), S. 342. 


4) F. Ringleb: Exakte Lösungen der Differentialgleichungen einer adiabatischen Gasströmung. Z. angew. Math. 
Mech. Bd. 20 (1940), S. 185 bis 198. 
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die Molenbroek-Tschapliginsche Transformation°) der gasdynamischen Differentialgleichun- 
gen. Kurz gesagt, wird danach die Stromfunktion y nicht als Funktion cartesischer Orts- 
koordinaten x, y, sondern als Funktion polarer Geschwindigkeitskoordinaten betrachtet, als 
welche Geschwindigkeitsbetrag q und Richtungswinkel # der Strömung gegen die =-Achse 
genommen werden. Die herauskommende lineare Differentialgleichung ist für die Herstellung 
von Beipielen besonders bequem. F. Ringleb fand bei seinen Strömungen die singulären 
Eigenschaften meiner Grenzlinien — bis auf die letzte®) — wieder. In einer anschließenden 
Arbeit ’) hat F. Ringleb nun die Anfänge einer allgemeinen Theorie der Grenzlinien entwickelt, 
wobei ihm wiederum die Molenbroek-Tschapliginsche Transformation sehr zustatten 
kam. Die Auseinandersetzung mit diesen einschlägigen Untersuchungen F. Ringleb's wird 
in den Abschnitten 4 und 5 der vorliegenden Abhandlung ihren Platz finden, wo das Verhalten 
der Stromlinien, des Druckgefälles und der Machschen Wellen an der Grenzlinie erörtert wird. 

Die erste Veranlassung zu meiner gegenwärtigen Untersuchung gab die Bemerkung, 
daß die Gleichung der Grenzlinie unmittelbar aus der Molenbroek-Tschapliginschen Trans- 
formation folgt. Die Gleichung der Grenzlinie bedeutet nämlich niehtsanderes als das Verschwinden 
einer Funktionaldeterminante, welche für die eineindeutige Abbildung eines @-y-Bereiches oder 
Ortsplanes auf einen g#-Bereich oder Geschwindigkeitsplan entscheidend ist. 
Gelegentlich der Untersuchung des Zusammenhanges von Orts- und Geschwindigkeitsplan 
wurden dann diejenigen Strömungen festgestellt, welche bei Anwendung der Molenbroek- 
Tschapliginschen Transformation der Betrachtung verlorengehen. Ferner sollte der Beweis 
geführt werden, daß die Grenzlinie allgemein Enveloppe einer Schar Machscher Wellen ist. 
Dieser Beweis war in der Ringlebschen Arbeit trotz eines dahingehenden Anspruches tat- 
sächlich nicht erbracht. Schließlich und vor allem wird bewiesen werden, daß irgendeine 
stationäre ebene Fortsetzuug der Strömung über eine Grenzlinie hinaus nicht möglich ist. 
Aus dieser Nichtfortsetzbarkeit folgt, daß die errechnete Potentialströmung sich bis zur Grenz- 
linie hin nicht verwirklichen läßt. Das theoretische Auftreten von Grenzlinien bedeutet des- 
halb einen Zusammenbruch der berechneten adiabatischen Potentialströmung schon strom- 
aufwärts vor der Grenzlinie. 


1. Die Molenbroek-Tschapliginsche Transformation. 


Zur Vereinfachung der Schreibweise verwenden wir im folgenden nur dimensionslose 

Geschwindigkeiten, die auf die größtmögliche Geschwindigkeit V, bezogen sind, welche 
F 5) 

beim Ausströmen ins Vakuum eintritt und gleich og ‚En ist, wo p, der Kesseldruck, 
o, die Kesseldichte und x das Verhältnis der spezifischen Wärmen ist. Mit « bezeichnen 
wir die dimensionslose Geschwindigkeitskomponente in der «-Richtung, mit ® die in der dazu 
senkrechten y-Richtung, mit q den Geschwindigkeitsbetrag, mit # den Winkel der Geschwindig- 
keit gegen die »-Richtung. Die Dichte denken wir uns durch Bezug auf die Kesseldichte 
dimensionslos gemacht. 


Dann gilt für die dimensionslose Dichte: 


1 


Die Schallgeschwindigkeit a ist bestimmt durch: 


a? 5 ( 


Wir merken gleich als Folgerung an: 


(a). 

5) P. Molenbroek: Über einige Bewegungen eines Gases bei Annahme eines Geschwindigkeit-potentials. Arch. 
Math. Phys. 1890, 2. Reihe, 9. Teil, S. 157 bis 195. A. Tschapligin: Über Gasstrahlen, Wissenschaftl. Ann. d. Univ. 
Moskau, Phys. Math. Kl. Bd. 21 (194), S. 1 bis 121 (russisch) wieder abgedruckt in der vollständigen Ausgabe seiner 
Werke. Über die Tsehapliginschen Untersuchungen findet sich ein Bericht in der Abhandlung von B. Demtehenko: 
Sur la relation entre la dynamique des fluides eompressibles et celle des fluides ineompressibles. Publie. Math, de 
l’Univers. de Belgrade Bd. 2 (1933), S. 85 bis 105. 


6) Vgl. meine Zuschrift an den Herausgeber in Z. angew. Math. Mech. Bd. 20 (1940), S. 295 bis 296. 


?) F. Ringleb: Über die Differentialgleichungen einer adiabatischen Gasströmung und den Strömungsstoß. 
Deutsche Mathematik Bd. 5 (1940), S, 377 bis 384. 
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Die beiden Grundgleichungen der stationären ebenen Potentialströmung sind die Gleichung 


der Wirbelfreiheit: 
dv du 
und die Kontinuitätsgleichung: 
dlor) 


Unter Einführung eines Potentials 9 und einer Stromfunktion y kann man zur Befriedigung 
der beiden Grundgleichungen setzen: 


dy . (ba), 

. (6b). 


Aus diesen beiden Differentialgleichungen folgt mit Benutzung von (3) die Differential- 
gleichung für allein: 


Diese Gleichung wollen wir fortan schlechthin als „Strömungsgleichung* bezeichnen. 
Beiläufig sei auch die Differentialgleichung für das Potential 9 angeführt: 


Die charakteristischen Grundkurven beider Differentialgleichungen genügen der Beziehung 
+2urdedy+la 
oder 
Diese charakteristischen Grundkurven sind im Überschallgebiet 4 > a reell; es sind die beiden 


Scharen Machscher Wellen, welche mit den Stromlinien den Machschen Winkel a=are sin r 


bilden. 
Potential- und Stromfunktion sind nach (6) definiert durch: 


Nach der Idee der Hodographenmethode sollen 9 und y nicht mehr als Funktionen der 
cartesischen Ortskoordinaten x, y, sondern der polaren Geschwindigkeitskoordinaten q, 9 auf- 
gefaßt werden. Um die Beziehung zwischen den alten unabhängigen Veränderlichen «, y 
und den neuen q, ® aufzustellen, löst man das System (10) nach d« und dy auf: 


f 


Indem @ und y jetzt als Funktionen von q und ® aufgefaßt werden, bekommt man unter 
Berücksichtigung von u=gc0sd, v»=gsin 


1 


Da x und y durch diese Beziehungen definiert sind, müssen nach dem Hauptsatz über 
Kurvenintegrale die Integrabilitätsbedingungen gelten: 


1 


1 2 
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Dies Differentialgleichungssystem entspricht dem Cauchy-Riemannschen bei inkom- 
pressiblen Strömungen (a>x), wobei man sich zu erinnern hat, daß dann In q und # isome- 
trische Koordinaten sind. Die beiden Differentialgleichungssysteme (12) und (13) bilden die 
Grundlage der Hodographenmethode der Gasdynamik. 


Die Beziehungen (12) zwischen Orts- und Geschwindigkeitskoordinaten formen wir für 
unsere Anwendungen mit Hilfe von (13) so um, daß 9 nicht mehr auftritt: 


sind . (14a), 


2 
cos + sind yy 


a? 
dy=,, sin . (14b), 


Aus dem Differentialgleichungssystem (13) kann man eine Differentialgleichung für y 
allein finden: 


Diese Differentialgleichung wollen wir nach ihrem bisher erfolgreichsten Bearbeiter „Tschapli- 
ginscheDifferentialgleichung“ nennen. Im Unterschied zur ursprünglichen Strömungs- 
gleichung ist die Tschapliginsche Differentialgleichung linear. 


Tschapligin hat eine Reihe von Partikularlösungen seiner Differentialgleichung ange- 
geben, die folgende Gestalt haben: 


Die Funktionen f, genügen einer hypergeometrischen Differentialgleichung. Für manche Fälle 
hat Tschapligin sogar die Konvergenz der Entwicklung einer durch Randbedingungen fest- 
gelegten Lösung nach solchen Partikularlösungen bewiesen. F. Ringleb hat eine Partikular- 


lösung dieser Art, nämlich ve ‚ erörtert und in dieser Strömung Grenzlinien festgestellt. 


Die charakteristischen Grundkurven der Tschapliginschen Differentialgleichung wollen 
wir als „charakteristische Hodographen“ bezeichnen. Sie gehorchen der Differential- 
gleichung 


Den Machschen Wellen als den charakteristischen Grundkurven im Ortsplan sind die 
charakteristischen Hodographen im Geschwindigkeitsplan zugeordnet. Für diese Zuordnung 
bedarf man nicht der Tschapliginschen Differentialgleichung, sie folgt vielmehr mittels 
der Charakteristikentheorie gleich aus der Strömungsgleichung. Die charakteristischen Hodo- 
graphen sind gut bekannt und werden in den sogenannten Charakteristikendiagrammen bei 
dem Prandtl-Busemannschen Verfahren‘) zur graphischen Ermittlung von Überschall- 
strömungen benutzt. 


Wir setzen das Integral 


Dann sind die beiden Scharen der charakteristischen Hodographen dureh: 
und 
bestimmt. Für die spätere Verwendung merken wir uns die Normalform der Tscehapli- 


ginschen Differentialgleichung für das Überschallgebiet an, die man bei Benutzung von & 
und n als unabhängigen Veränderlichen erhält: 


Die Koeffizienten dieser Differentialgleichung, die noch als Funktionen von g geschrieben sind, 
hat man sich natürlich durch & und n ausgedrückt zu denken. 


s) Vgl. Handbuch der Experimentalphysik Bd. 4, Teil 1, S, 424, Leipzig 1931. 
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2. Durch die Molenbroek-Tschapliginsche Transformation verlorene Strömungen. 


Der Übergang von der Strömungsgleichung zur Tschapliginschen Differentialgleichung 
kommt auf die Benutzung der Geschwindigkeitskoordinaten 4, ® anstatt der cartesischen Orts- 
koordinaten x, y als unabhängiger Veränderlicher für die Stromfunktion hinaus. Also müssen 
sich die Ortskoordinaten &, y als eindeutige Funktionen von q und d ausdrücken lassen. 
Wenn dies nicht der Fall ist, so ist die Transformation nicht gestattet. 

Die eindeutige Umkehrbarkeit der Funktionen q (x, y) und # (x, y) ist nieht mehr ge- 
währleistet, wenn die Funktionaldeterminante 
determinante längs einer Kurve verschwände, so könnten in deren Umgebung zu demselben 
q-d-Bereich verschiedene x-y-Bereiche gehören. Dann wäre nur eine entsprechende Abgren- 
zung der «-y-Bereiche nötig. Dagegen ist bei „identischem“ Verschwinden der Funktional- 
determinante in einem «-y-Gebiete mit von 0 verschiedenem Flächeninhalt die Molenbroek- 
Tsehapliginsche Transformation überhaupt nicht mehr anwendbar. Wenn die Strömungs- 
gleichung in einem solchen Falle Integrale hat, so gehen diese beim Übergang zur Tschapli- 
ginschen Differentialgleichung verloren. 

Zur Aufstellung dieser verlorenen Integrale benutzen wir die Kontinuitätsgleichung (5) 
unter Einführung des Ausdruckes (1) für o und die Gleichung der Wirbelfreiheit (4): 


verschwindet. Wenn diese Funktional- 


ou, du or 
du 
9) 


Bei identischem Verschwinden von >a,n) wird # eine Funktion von q allein: 9 (q). In den 


beiden Gl. (20) setzen wir u—=gcosd, v—qsind und erhalten dann: 


Ü 
— cos — ga?sin + g})sind + cos 2 =0. . .(2la), 


Damit dies homogene lineare Gleichungssystem für 24 und r nicht triviale Lösungen 


hat, d. h. damit q nicht identisch konstant wird, muß nach der Determinantenbedingung 
-(3-2)=0 
sein. Diese Beziehung für die gesuchte Funktion 9 (g) liefert zusammen mit (21b) den 
Neigungstangens für die Linien konstanter Geschwindigkeit g: 


dd 
sin d 4 cos 
Da nach (22) 
dd 
=+tciga 


sein soll, wo a der Machsche Winkel ist, so erhält man schließlich für den erwähnten 
Neigungstangens tg (9 +.a), d. h. die Linien mit konstantem q fallen mit einer Schar 
Machscher Wellen zusammen. Da d und a als Funktionen von g allein auf einer solchen 
Schar ebenfalls konstant sind, ist der Neigungstangens jeder Machschen Welle dieser 
Schar konstant. Diese Machschen Wellen sind also Geraden. y ist natürlich längs dieser 
Geraden veränderlich, während 4 und # konstant sind, so daß y unmöglich als Funktion von q 
und # betrachtet werden kann, wie die Molenbroek-Tschapliginsche Transformation 
es will. 
Diese „verlorenen Strömungen“ sind sehr geläufig, sie treten als Überschall- 
strömungen bei einseitiger Begrenzung auf. Daß die Schar geradliniger Machscher Wellen 
einer solchen Überschallströmung auch eine Enveloppe bilden kann, erkennt man bei Annahme 


en 


>) 
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einer Begrenzung durch eine Wand, die zunächst eben, dann mit der Höhlung nach der 
Strömung zu gekrümmt ist. Eine Enveloppe Machscher Wellen wird in einem solchen Falle 
bei Verzögerung und Kompression gebildet, während in den bekannten Beispielen zu der 
folgenden allgemeinen Theorie (vgl. Abschnitt 5) die Enveloppenbildung bei Beschleunigung 
und Expansion auftritt (s. Bild 1 und 2). 


Bild 1. Ankommende Strömung mit Grenz- 
linie nach W. Tollmien, Z. angew. & 
Math. Mech. Bd. 17 (1937), S.131. Die dieken 
Linien bezeichnen Stromlinien, die dün- 


Bild2. Druck p über der untersten 
Stromlinie der Strömung von 
Bild 1senkrecht aufgetragen.p, ist 
der Kesseldruck. Durch Striche- 


nen Machsche Wellen, zu denen auch die 4 lung ist die Stelle hervorgeho- 
Grenzlinie ganz rechts gehört. Auf der ben, wo Schallgeschwindigkeit 
kurzgestrichelten Linie ist die Strömungsgeschwindigkeit herrseht. 


gleich der Schallgeschwindigkeit. = 1,405. 


3. Mehrdeutigkeit von Strömungen an der Grenzlinie. 


Fortan beschäftigen wir uns nur noch mit Strömungen, die wir uns aus regulären Lösungen 
y(q,9) der Tschapliginschen Differentialgleichung gewonnen denken können. Die 
Stromfunktion y und die Ortskoordinaten «, y seien als eindeutige Funktionen von q und # 
gegeben. Es erhebt sich nun die Frage, ob man diese Funktionen x (q, 9) und 9 (g, ®) immer 
eindeutig umkehren kann. Diese Frage ist das Gegenstück zu der im vorigen Abschnitt. Wir 
wollen jetzt wissen, wann wir mit der Möglichkeit rechnen können, daß zu einem Werte- 
paar x, y verschiedene Wertepaare 4, 9 gehören, so daß eine Mehrdeutigkeit der Strömung 
im Ortsplan (x-y-Bereich) statthat. 

Die Eindeutigkeit der Umkehrung ist nicht mehr gewährleistet, wenn die Funktional- 


determinante verschwindet. Diese Funktionaldeterminante wird nach (14) berechnet: 


Zunächst zeigen wir, daß diese Funktionaldeterminante nicht in einem q-d-Bereich mit 
von O0 verschiedenem Flächeninhalt verschwinden kann. Wäre es an dem, so würde die 
Tschapliginsche Differentialgleichung überschüssige Integrale liefern, welchen nicht Lösungen 
der Strömungsgleichung entsprechen. 

Zum Beweise nehmen wir identisches Verschwinden der Funktionaldeterminante an. 
In dem betreffenden Gebiete würde dann die Differentialgleichung 


ya 


gelten, was höchstens im Überschallgebiet möglich wäre. Führen wir als unabhängige Ver- 
änderliche & und „7 nach (18) ein, so bekommt man daraus die Differentialgleichung: 


mit den beiden Möglichkeiten y:—0 oder yn=0. Beide Male erhält man einen Widerspruch 
zu der Tschapliginschen Differentialgleichung in ihrer Normalform (19). 

Im Unterschallgebiet 9<a kann die Funktionaldeterminante nur in einem Punkte 
yya=0, ys—=0 verschwinden. Daß ein solcher Punkt ein Staupunkt ist, liegt nahe und wird 
auch von Ringleb behauptet. Sein Beweis ist aber nicht stichhaltig, da er in diesem Falle 
noch das Nichtverschwinden der Funktionaldeterminante ar voraussetzt. Wir schließen 
so: Das Differential 

d d 4 
verschwindet im Punkte y,=0, y,=0. Andererseits ist dies Differential nach (10b): 


dy=o(udy—-vde). 
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Da dy zu d« im beliebigen Verhältnis stehen kann, müssen «=0, v—=0 sein, damit dy 
immer verschwindet. Also ist der Punkt y,=0, ys=0 ein Staupunkt. Das gleichzeitige 
Verschwinden von y, und y'; ist hiermit vollständig erledigt und kann im folgenden stets 
ausgeschlossen werden. 

Im Überschallgebiet q9>a besteht noch die Möglichkeit, daß die Funktional- 
determinante längs einer Kurve des Geschwindigkeitsplanes verschwindet, wofür auch Bei- 
spiele vorliegen. Eine solche Kurve wollen wir „Grenzhodograph“ nennen, während die 
zugehörige Bildkurve im Ortsplan eine Grenzlinie ist. Es ist nämlich wegen des Ver- 
schwindens der Funktionaldeterminante damit zu rechnen, daß in der Umgebung einer so defi- 
nierten Grenzlinie zu den x-y-Werten verschiedene 4-9Werte und damit auch verschiedene 
y-Werte gehören. Die Integralfläche y (x, y) würde dann an der Grenzlinie eine Rückkehr- 
kante haben. Um «diese Vermutung zu bestätigen, werden wir das Verhalten der Strom- 
linien y= const, d. i. der Höhenlinien der Integralfläche y (x, y) erörtern. 

Eine maßgebende Funktionaldeterminante der Molenbroek - Tschapliginschen 

_ Transformation hat uns also zwanglos auf die Gleichung des Grenzhodographen: 


geführt. Da die weiteren Untersuchungen zu einem guten Teil in dem bequemen Geschwin- 
digkeitsplan angestellt werden, wollen wir die Gl. (24) für das Geschwindigkeitsbild der 
Grenzlinie als „Garenzgleichung“* schlechthin bezeichnen. Im folgenden werden wir aber 
diese Grenzgleichung nicht voraussetzen, sondern sie jeweils aufs neue aus den eingangs er- 
wähnten, von Beispielen her bekannten singulären Eigenschaften ableiten, welche die Strom- 
linien, das Druckgefälle und die Machschen Wellen an der Grenzlinie aufweisen. 


4. Stromlinien und Druckgefälle an der Grenzlinie. 
Längs einer Stromlinie y = const gilt: 


also 


Setzen wir dies in (14) ein, so bekommen wir für die Komponenten des Linienelementes 
einer Stromlinie: 


sin 1 


q ist also hier als Parameter längs der Stromlinie genommen. Aus (26) ergibt sich das 
Linienelement einer Stromlinie: 


1 1 1 
Für das Druckgefälle 


in Stromrichtung erhält man aus der letzten Gl. (27): 


ya— ® 


Damit das Druckgefälle unendlich wird, muß wieder die Grenzgleichung (24) gelten. Daß 
nicht etwa ein Verschwinden von y; (im Zähler von (28)) diesen Schluß umstößt, werden 
wir bald durch eine kleine Diskussion der Grenzgleichung sehen. 
®) Wenn ein Verschwinden von wg befürchtet wird, so kann ınan die Beziehung dy—-- "Das benutzen; ein 
a 


gleichzeitiges Verschwinden von wg und Ygq kommt ja nach einer früheren Bemerkung nieht mehr in Frage. Man 
würde auch mit dieser Beziehung bei den nachfolgenden Betrachtungen genau so gut wie mit (25) auf die Grenz- 
gleichung kommen, 
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Als Bedingung für das Auftreten singulärer Punkte FR -5) der Stromlinien erhalten 


wir nach (26) ebenfalls die Grenzgleichung. 

Die bisherigen Ergebnisse dieses Abschnittes finden sich bereits bei F.Ringleb, wenn 
auch mit etwas anderer Herleitung. 

Wir bemerken an der Grenzgleichung, daß y, für die Grenzlinie nur dann 0 werden 
kann, wenn q=a ist; das gleichzeitige Verschwinden von y, und y; wurde ja ausgeschlossen. 
y; kann für die Grenzlinie nur verschwinden, wenn a=0 ist, was uns wegen des dann be- 
stehenden Vakuums auch nichts angeht. 

Die Art der singulären Punkte der Stromlinien auf der Grenzlinie erkennt man leicht, 
wenn man sich daran erinnert, daß die Richtung der Stromlinien durch den Wert der 
Molenbroek-Tschapliginschen unabhängigen Veränderlichen 9 in dem betreffenden 
Punkte der Grenzlinie eindeutig gegeben ist. Die Stromlinien könnten also entweder knick- 
frei die Grenzlinie durchlaufen oder sie kehren mit derselben Richtung unter Bildung einer 
Schnabelspitze um. Im allgemeinen liegt der zweite Fall vor, da der Ausdruck 


1 1 2 


an der Grenzlinie in der Regel eine einfache Nullstelle hat, während ys nach der obigen Be- 
merkung dort nicht verschwinden kann. Dann erleiden dx und dy nach (26) an der Grenz- 
linie einen Vorzeichenwechsel, die Stromlinien kehren um. 


Daß der Ausdruck (29) im allgemeinen eine ein- 
fache Nullstelle hat, kann man an seiner Ableitung 
längs einer Stromlinie nach einem Parameter, als welchen 
wir q wählen, erkennen. Unter Berücksichtigung von (25) 
erhält man für diese Ableitung 


1 1 


44 


Die Tsehapliginsche Differentialgleichung und die 
Grenzgleichung erzwingen nicht unmittelbar das Ver- 
Bild 3. Ankommende Strömung wie in Schwinden dieser Ableitung (30), so daß der Ausdruck (29) 
Bild 1 mit Rückkehrströmung. Die Strom- jm allgemeinen sein Vorzeichen an der Grenzlinie wech- 
linien der ankommenden Strömung sind It Die bek t Beispiel . d .] . 
ausgezogen, die der Rückkehrströmung lang. Sell. ıe bekann en e eispiele zeigen denn auch eine 
gestrichelt. Auf der kurzgestrichelten Linie Rückkehr der Stromlinie. 


herrscht für die ankommende Strömung 
Schallgeschwindigkeit. Die an der Grenzlinie „ankommende Strömung“ 


würde also formal durch die „Rückkehrströmung* 
fortgesetzt (s. Bild 3). Einen physikalischen Sinn hat diese Art der Fortsetzung natür- 
lich nicht. 


5. Machsche Wellen an der Grenzlinie. 

Von den Beziehungen der Grenzlinie zu den Machschen Wellen weisen wir zunächst 
die Tatsache nach, daß die Grenzlinie selbst eine Machsche Welle ist. 

Zu dem Zweck transformieren wir die Differentialgleichung (9) der Machschen Wellen 
auf die Veränderlichen y, q, indem wir u=gc0osd, v=qsin® setzen und und dy 
aus (14) entnehmen. Diese Rechnung ergibt unmittelbar 

a? 1 a? 
Nach einer kleinen Umordnung erhält man die transformierte Differentialgleichung der 
Machschen Wellen in der endgültigen Form: 
1 
v3 


1 1 
Man bekommt beim Verschwinden der ersten Klammer von (31) diejenigen Machschen 


Wellen, die im Geschwindigkeitsplan durch die charakteristischen Hodographen (16) dar- 
10° 
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gestellt sind. Außerdem erhält man aber noch als Machsche Welle die Grenzlinie, deren 
Gleichung im Geschwindigkeitsplan 


ist. 

An den Beispielen war seinerzeit von mir festgestellt worden, daß die Grenzlinie eine 
Schar Machscher Wellen einhüllt. Damit ist keineswegs gesagt, daß der Grenzhodograph 
als Geschwindigkeitsbild der Grenzlinie die Bilder einer Schar Machscher Wellen im Ge- 
schwindigkeitsplan, also eine Schar charakteristischer Hodographen einhüllt. Tatsächlich 
kennt man ja die charakteristischen Hodographen sowohl formelmäßig (18) als auch anschau- 
lich aus den sogenannten Charakteristikendiagrammen des Prandtl-Busemannschen Ver- 
fahrens, wo q und # als Polarkoordinaten der charakteristischen Hodographen aufgetragen 
sind. Die charakteristischen Hodographen besitzen danach keine Enveloppe außer bei der 
Ausströmungsgeschwindigkeit ins Vakuum (q=1). 

Während der eingangs durchgeführte Beweis, daß die Grenzlinie eine Machsche Welle 
ist, sich genau so bei F. Ringleb findet, ist der zum Schluß seiner Arbeit unternommene 
Beweisversuch, daß die Grenzlinie eine Schar Machscher Wellen einhüllt, mißglückt, da er 
im Geschwindigkeitsplan bleibt. Ringleb kombiniert mit der Grenzgleichung die Beziehung 


die auf einer Stromlinie gelten würde. Dann würde die Grenzgleichung in die Differential- 
gleichung der charakteristischen Hodographen (16) übergehen. Daraufhin wird der ver- 
blüffende Schluß gezogen, daß die Grenzlinie Einhüllende der Machschen Wellen und der 
Stromlinien ist. Tatsächlich kann die Grenzlinie als Machsche Welle nicht mit einer Strom- 
linie zusammenfallen, es sei denn im Vakuum, wo dann auch wirklich die charakteristischen 
Hodographen eine altbekannte Enveloppe besitzen. 

Das Verhalten einer Schar Machscher Wellen an der Grenzlinie untersuchen wir, indem 
wir die Linienelemente der Machschen Wellen aufstellen, also in den Ortsplan zurück- 
gehen. Da für eine Schar Machscher Wellen im Geschwindigkeitsplan nach (16) 


gilt, worin beide Wurzelvorzeichen möglich sind, je nachdem welche Schar ausgewählt ist, 
so erhält man nach (14) als Komponenten der Linienelemente einer Schar Machscher Wellen 


cos ya) +[eosdgy, — sinds] 

dy= — sin eos +Hsin day + cos ws] dgq. 


Nach einer kleinen Umordnung bekommt man: 


1 


1 
dy= 


Die beiden vorderen Klammern in (33a) und (33b) könnten gleichzeitig nur bei “ —=0 ver- 
schwinden, was bei Vorhandensein Machscher Wellen sinnlos ist. Als Ort der singulären 
Punkte (4 = 5) einer Schar Machscher Wellen ergibt sich die Grenzlinie, deren Gleichung 
im Geschwindigkeitsplan jetzt in der Form: 


i 


auftritt. Hiermit ist wiederum eine neue Definition der Grenzlinie gewonnen. 


Da die Wurzel in der Gl. (34) auf einer bestimmten Grenzlinie nur ein ganz bestimmtes 
Vorzeichen haben kann, so geht aus dem Beweisgang hervor, daß jeweils nur eine Schar 
Machscher Wellen auf der Grenzlinie singuläre Punkte haben kann, die wegen der bekannten 
Eigenschaften der Stromlinien an der Grenzlinie nur Schnabelspitzen sein können. 
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Für den Beweis, daß die andere Schar Machscher Wellen von der Grenzlinie ein- 
gehüllt wird, benötigen wir die grundlegende Tatsache, daß das Geschwindigkeitsbild der 
Grenzlinie kein charakteristischer Hodograph ist. 

Würde der Grenzhodograph ein charakteristischer Hodograph sein, so müßten für ihn 
gleichzeitig die beiden Gleichungen 


gelten. Durch Elimination von (=-5) würde sich daraus als notwendige Bedingung für 
den Richtungstangens'”) des (renzhodographen die Beziehung 


oder 
d# 


ergeben. Nun ist — Ya der Richtungstangens eines Stromlinienbildes. Also würde die erste 


Möglichkeit (35a) Zusammenfallen der Grenzlinie mit einer Stromlinie bedeuten; denn die 

Richtungsübereinstimmung überträgt sich aus dem Geschwindigkeitsplan in den Ortsplan, ob- 

wohl das Umgekehrte in der Nähe der Grenzlinie nicht mehr gilt. Die Grenzlinie als Mach- 

sche Welle kann aber mit einer Stromlinie nur im Vakuum zusammenfallen. Bei der zweiten 


Richtungstangens eines charakteristischen Hodographen sein würde, falls das von r gelten 


Möglichkeit (35b) gehen wir von der nach (16) folgenden Bemerkung aus, daß auch 


würde, was wir bis zur Aufweisung eines Widerspruchs gerade annalımen. Nach (35b) müßte 


also jetzt der Richtungstangens % eines charakteristischen Hodographen mit dem Rich- 


Y@ eines Stromlinienbildes zusammenfallen. Nach der Übertragbarkeit einer 
+ 


tungstangens — 


Richtungsübereinstimmung aus dem Geschwindigkeitsplan in den Ortsplan hätte man auf 
Übereinstimmung der Richtung von Machschen Wellen und Stromlinien zu schließen, was 
nur im Vakuum möglich wäre. Der innere Grund für die Notwendigkeit, die beiden Mög- 
lichkeiten (35a) und (35b) zu berücksichtigen, liegt darin, daß es zwei Scharen Machscher 
Wellen in der Umgebung der Grenzlinie gibt. 

Die zur Diskussion stehende lückenlos den Geschwindigkeitsplan überdeckende Schar 
charakteristischer Hodographen wird also durch den Grenzhodographen geschnitten. Über- 
tragen wir nun diese charakteristischen Hodographen aus dem Geschwindigkeitsplan als 
Machsche Wellen in den Ortsplan, so treffen diese die Grenzlinie. Da die Grenzlinie als 
Machsche Welle genau so gut den Machschen Winkel mit den Stromlinien einschließt wie 
die betreffende Schar Machscher Wellen, so hüllt die Grenzlimie diese Schar Machscher 
Wellen ein. Dabei ist, wie wir schon früher bemerkten, von den beiden Scharen Machscher 
Wellen diejenige genommen, welche keine singulären Punkte an der Grenzlinie hat. 

Die Anwendung der Charakteristikentheorie auf die Strömungsgleichung (7) 
liefert im Geschwindigkeitsplan nur die charakteristischen Hodographen, dagegen kein Ge- 
schwindigkeitsbild der Grenzlinie, obwohl die Grenzlinie als Machsche Welle eine charakte- 
ristische Grundkurve ist. Dies ist jetzt nicht weiter verwunderlich, da die Grenzlinie als 
Enveloppe Machscher Wellen keine neuen charakteristischen Linienelemente auf der Integral- 
fläche y (x, y) einführt, die nicht schon den charakteristischen Hodographen zugeordnet sind, 
sondern nur eine andere Zusammenfassung dieser charakteristischen Linienelemente vornimmt. 

Beiläufig wollen wir bemerken, daß man aus der Differentialgleicehung (9) der Mach- 
schen Wellen auch mit der Legendreschen Transformation eine kennzeichnende Beziehung 
für die Grenzlinie erhalten kann. Diese Legendresche Transformation linearisiert die 
Differentialgleichung (8) für das Potential 9. Als neue unabhängige Veränderliche werden 


s 10) Wir denken uns jetzt einmal der einfacheren Ausdrucksweise wegen q und # als cartesische Koordinaten 
aufgetragen. 
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eingeführt, als neue abhängige Veränderliche „=ur+ry-—g. Die Um- 
rechnung der Differentialrechnung (9) der Machschen Wellen auf die neuen Veränderlichen 
ergibt 


la? — dw — 2urdude+(a® — —ZunXer]=0. 


Während das Verschwinden der ersten Klammer dieser Gleichung die charakteristischen 
Hodographen definiert, kennzeichnet 


Zus — 


die Grenzlinie. Dieselbe Beziehung würde aus dem Verschwinden der Funktionaldeterminante 
(u, 
wir aber nieht unternommen. 


folgen. Einen Aufbau der Grenzlinientheorie von dieser neuen Grundlage aus haben 


6. Nichtfortsetzbarkeit von Strömungen mit Grenzlinie. 


Bis zum Beweis des Gegenteils nehmen wir vorläufig an, daß die berechnete ankommende 
Strömung bis zur Grenzlinie verwirklicht werden kann. Wir fragen nach der physikalischen 
Fortsetzung dieser Strömung. Wie in der ganzen Abhandlung beschränken wir uns stets auf 
stationäre ebene Strömungen. 

Zwei Möglichkeiten der Fortsetzung haben wir zu unterscheiden. Die „Fortsetzungs- 
strömung“ kann entweder ohne Unstetigkeit in den Geschwindigkeiten an die ankommende 
Strömung anschließen oder die Grenzlinie wird Träger von Geschwindigkeitssprüngen. 
Geschwindigkeitssprünge bedeuten Unstetigkeiten in den ersten Ableitungen von y' nach 
den Ortskoordinaten. 

Wir betrachten zuerst die zweite Möglichkeit. In dem Falle würde die Grenzlinie Front 
eines Verdichtungsstoßes sein; denn eine Wirbelschicht können die Geschwindigkeitssprünge 
nicht bilden, da die Grenzlinie keine Stromlinie ist. Nach der Verdichtungsstoßtheorie 
(vgl. etwa die Busemannsche Stoßpolare) ist bei gegebener Anströmgeschwindigkeit der 
Geschwindigkeitssprung eindeutige Funktion des Winkels der Stoßfront mit der An- 
strömungsrichtung. Die Grenzlinie als Stoßfront würde mit der Anströmungsrichtung den 
Machschen Winkel bilden. Bei diesem Winkel verschwindet aber nach der Verdichtungs- 
stoßtheorie der Geschwindigkeitssprung. 

Bei der übrigbleibenden Möglichkeit sind also beim Durchgang durch die Grenzlinie 
nicht nur die Stromfunktion y, sondern auch die Geschwindigkeitsgrößen q und # sowie 
Druck und Dichte stetig. Die Strömung müßte beim Überschreiten der Grenzlinie adiabatisch 
und darum wirbelfrei bleiben. Man hat dann noch zwei Fälle zu unterscheiden, erstens 
Fortsetzung durch verlorene Strömungen (vgl. Abschnitt 2), zweitens Fortsetzung durch 
Strömungen, die aus der Tschapliginschen Differentialgleichung gewonnen werden können. 

Um den ersten Fall zu behandeln, erinnern wir uns, daß bei den verlorenen Strömungen 
die Kurven mit konstantem 4 Geraden sind, auf denen gleichzeitig 9 konstant ist. Ist nun 
auf der Grenzlinie oder jedenfals auf einem Stück der Grenzlinie 4 konstant, so müßte dort 
wegen des erforderlichen stetigen Anschlusses der ankommenden Strömung an eine verlorene 
Strömung auch ® konstant sein und die Grenzlinie müßte dort geradlinig verlaufen, was in 
jedem Einzelfalle leicht zu widerlegen wäre. Wir können aber auch allgemein schließen. 
Da die ankommende Strömung nämlich bis zur Grenzlinie einschließlich voraussetzungsgemäß 
aus der Tschapliginschen Differentialgleichung gewonnen werden kann, so können auf 
der Grenzlinie q und Ö nicht gleichzeitig konstant sein, weil y dann nicht als Funktion von 
q und Ö aufgefaßt werden könnte. Wenn nun andererseits q auf der Grenzlinie nicht konstant 
ist, so muß wegen des stetigen Anschlusses # als Funktion von q wiederum aus der Theorie 
der verlorenen Strömungen, also jetzt aus Gl. (22), bestimmt sein. Diese Gleichung ist aber 
nichts anderes als die charakteristische Hodographengleichung (16), die für die Grenzlinie 
gerade nicht gilt, wie wir im Abschnitt 5 bewiesen. 

Als letzter Ausweg ist demnach nur noch die Fortsetzung durch eine Strömung ge- 
blieben, die aus einer Tschapliginschen Differentialgleichung gewonnen werden kann, Wir 
suchen jetzt die Fortsetzungsströmung aus einem Anfangswertproblem der Tschapligin- 
schen Differentialgleichung -zu bestimmen. Die Anfangswerte müssen natürlich auf dem 
Grenzhodographen angegeben werden. Auf allfällige Randeinflüsse kommen wir noch zu 
sprechen. 

Auf den Grenzhodographen führen wir einen Parameter o ein. Als Funktionen dieses 
Parameters sind y (0), 4(9), ®(0), 0o(6) aus der ankommenden Strömung wegen des stetigen 
Anschlusses von Stromfunktion, Geschwindigkeiten und Dichte auch für die Fortsetzungs- 
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strömung bekannt. Ferner kennen wir die Koordinaten der Grenzlinie = (0), y(o). Da nach 
(10b) in der ankommenden Strömung 


dy=oudy-ordzx 


sein muß, so besteht zwischen diesen sechs Funktionen die Beziehung 


„dx 


Als erste Anfangsbedingung ist y(o) auf dem Grenzhodographen gegeben. Da der Grenz- 
hodograph kein charakteristischer Hodograph ist, so müssen zur vollständigen Stellung des 
Anfangswertproblems noch die beiden Ableitungen y, und ys für die Fortsetzungsströmung 
auf dem Grenzhodographen gegeben sein. Däbei ist zu beachten, daß man noch nicht weiß, 
ob y, und y, beim Übergang von der ankommenden Strömung zur Fortsetzungsströmung 
stetig bleiben. Aus diesem Grunde darf die Gleichung 


ya ar 


welche für den Grenzhodographen der ankommenden Strömung gilt, nicht benutzt werden. 
Zur Verfügung stehen uns die beiden Gl. (14), welche für jede aus der Tschapliginschen 
Differentialgleichung erzielbare Strömung, also auch die vermutete Fortsetzungsströmung, 
gelten. Diese beiden Gleichungen liefern für den Anfang der Fortsetzungsströmung 


do 


cos (1 sind +[eosdgy, — sind ws] . (37a), 


Außerdem gilt noch die Gleichung 


dy dq 


Wegen der Bedingung (36) folgt diese Gl. (35) aus den beiden vorhergehenden, so daß (38) 
nicht weiter berücksichtigt zu werden braucht. Aus den beiden Gl. (37) lassen sich nun die 
gesuchten Anfangswerte y, und y; für die Fortsetzungsströmung eindeutig bestimmen. 
Die Determinante dieses Gleichungssystems für y, und ys wird nämlich 


da) 
(ao) 


Verschwände diese Determinante, so würde der Grenzhodograph ein charakteristischer Hodo- 
graph sein, was er bewiesenermaßen nicht ist. 

Mit demselben Recht wie für die Fortsetzungsströmung könnte man die obige Be- 
trachtung auch für die ankommende Strömung durchführen, und man erhielte die gleichen 
Werte y, und y,;, da die benutzten Funktionen y (0), 4.(0), (0), 0 (0), x (0), y(0) für beide 
Strömungen die gleichen sind. Die ankommende und die Fortsetzungsströmung stimmen also 
nicht nur in y, sondern auch in y, und y, längs des Grenzhodographen überein. 


Bild 4. Abhängigkeitsgebiete der ankommen- 
den Strömung von Bild 1 und der vermuteten 
Fortsetzungsströmung im Gesehwindigkeitsplan. 


T Die dieke Linie stellt den Grenzhodographen 
$ dar. Daß q auf dem Grenzhodographen dieses 
Bildes konstant ist, liegt nur an dem besonderen 


Beispiel. Die dünnen Linien sind eharakteristi- 
sche Hodographen. Auf der gestrichelten Linie 
ist 


[7280-#] 


Wir betrachten nun einen solchen Bereich des Parameters o, daß das entsprechende 
Stück der Grenzlinie x (0), y(0) ganz innerhalb der Umrandung liegt, welche für die Strömung 
vorgesehen ist. Über dem entsprechenden Stück des Grenzhodographen denken wir uns die 
beiden „Abhängigkeitsgebiete“ gemäß der Theorie linearer hyperbolischer Differential- 
gleichungen konstruiert (s. Bild 4). Das sind zwei krummlinige Dreiecke) über dem Grenz- 

11) In Bild 4 ist die Spitze eines dieser Dreiecke durch die gestrichelte Linie q=a weggeschnitten, da nach 


Unterschreitung der Schallgeschwindigkeit die Tschapliginsche Differentialgleichung elliptisch ist. Die weggeschnittene 
Spitze gehört nieht mehr zum Abhängigkeitsgebiet im Sinne der Theorie hyperbolischer Differentialgleichungen. 
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hodographen als gemeinsamer Grundseite, deren andere Seiten durch charakteristische Hodo- 
graphen gebildet werden. In diesen Abhängigkeitsgebieten ist nach der Theorie linearer 
hyperbolischer Differentialgleichungen y (q, 9) vollständig durch die Anfangsbedingungen auf 
dem Grenzhodographen bestimmt. Einflüsse irgendwelcher Randbedingungen machen sich 
in diesen Abhängigkeitsgebieten noch nicht bemerkbar. 

Nun kennen wir bereits zwei Lösungen, welche sich mit den gleichen Werten y, y,, y's 
auf dem Grenzhodographen aneinanderschließen; sie entsprechen der ankommenden Strömung 
und der Rückkehrströmung, die ja aus ein und demselben Integral der Tschapliginschen 
Differentialgleichung gewonnen wurden. Wegen der eindeutigen Bestimmtheit der Lösung 
durch die Anfangswerte y, y,, ys auf dem Grenzhodographen muß die vermutete Fort- 
setzungsströmung in dem betreffenden Abhängigkeitsgebiet des Geschwindigkeitsplanes mit 
der Rückkehrströmung zusammenfallen. Im Ortsplan würde dies Übereinstimmung der Fort- 
setzungsströmung mit der Rückkehrströmung in der Nachbarschaft der Grenzlinie bedeuten, 
gleichgültig wie die Umrandung weitergeführt wird. 

Damit entfällt auch .die letzte noch offene Möglichkeit einer Fortsetzung der an der 
Grenzlinie ankommenden adiabatischen Potentialströmung durch eine stationäre ebene Strömung. 
Deshalb läßt sich die berechnete adiabatische Potentialströmung bis zur Grenzlinie hin nicht 
verwirklichen, sondern bricht schon stromaufwärts zusammen. Die physikalische Bedeutung 
dieses Sachverhaltes wird freilich erst dann endgültig klar werden, wenn es gelingt, die 
tatsächliche Strömung innerhalb einer Umrandung herzustellen, die bis zum Ort der Grenz- 
linie mit dem Rande der Grenzlinienströmung übereinstimmt. 280 


Unstetigkeitsstellen in einem bewegten Kontinuum. 
Von W. Kucharski in Berlin. 


In Entwickelung eines von Christoffel uw. a. angegebenen Verfahrens 
werden die Grundgleichungen aufgestellt, die an einer Unstetigkeitsstelle 
der mafigebenden physikalischen Größen in einem Kontinuum gelten. An 
einem Beispiel wird gezeigt, daß sich hieraus mit häufig zulässigen und 
üblichen Vernachlässigungen Beziehungen für die Sprunggrößen ergeben, 
in denen nur diese selbst und von vornherein gegebene Funktionen des 
betreffenden Problems vorkommen, nicht aber die Rand- oder Anfangs- 
bedingungen. Ein Hauptergebnis besteht darin, daß sich bei den erwähnten 
Vereinfachungen die einzelne Unstetigkeitsstelle stets mit der „Schallge- 
schwindigkeit des Problems“ fortpflanzt. 


1. Einleitung. In den folgenden Abschnitten soll an einem einfachen Beispiel gezeigt 
werden, wie man gewisse Aussagen über die maßgebenden physikalischen Größen gewinnen 
kann, wenn sich in einem bewegten Kontinuum eine Unstetigkeitsstelle z. B. für die Spannung, 
die Dehnung und die Geschwindigkeit befindet. Es handelt sich dabei um die Anwendung 
und Weiterbildung einer Methode, die, soweit ich aus der von mir benutzten Literatur ersehe, 
von Christoffel für lineare Saitenschwingungen und von Hugoniot für Vorgänge in 
strömenden Gasen benutzt wurde‘). Es wird sıch zeigen, daß man unter gewissen häufig 
erfüllten Voraussetzungen eine einfache Auskunft über die Geschichte der Unstetigkeitsstelle 
erhält, ohne daß die vollständige Lösung der betreffenden Differentialgleichung bekannt zu 
sein braucht. 

Die Problemstellung wurde durch Untersuchungen an stabförmigen Ventilkörpern an- 
geregt, deren Querschnittsverlauf in einer für die Weiterbehandlung des Problems brauchbaren 
Weise analytisch nur mit großen Umständlichkeiten faßbar war, so daß die vorstehend an- 
gedeutete Fragestellung, ob man nicht ohne Kenntnis der vollständigen Lösung Daten über 
die Größenordnung der auftretenden Spannungen und über ihre Beeinflussung durch den 
(uerschnittsverlauf erhalten kann, sich fast von selbst ergab. 

Mathematisch gesehen hängt diese Methode mit der Theorie der Charakteristiken der be- 
treffenden hyperbolischen partiellen Differentialgleichungen zusammen’), von denen, wie bekannt, 
mehrere Integralflächen ausgehen und auf denen auch Unstetigkeiten gewisser Art bestehen 
können. Man findet hierüber im 2. Band von Courant-Hilbert, Methoden der mathematischen 


1)S. Frank-Mises: Differential-Gleichungen der Physik: Wien 1935, Band II, S. 323 usw.; ferner: Webster- 
Szegö: Partielle Differentialgleichungen der mathematischen Physik, Leipzig und Berlin 1959. : 


2) Hierauf wies Herr Hamel bei einem Seminarvortrag in dankenswerter Weise hin. 
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Physik, die allgemeinen Grundlagen. Diese weitreichenden Untersuchungen wurden bei der 
folgenden Arbeit, die seit etwa 8 Jahren im Manuskript vorlag, nicht benutzt. Die Veröffent- 
liehung dürfte sich rechtfertigen, denn diese Methoden scheinen — jedenfalls in der tech- 
nischen Mechanik — wenig bekannt zu sein; auch scheinen derartige Probleme in der letzten 
Zeit häufiger aufzutreten. Ich lasse es daher auch bei der ursprünglichen, eng an die physi- 
kalische Anschauung angelehnten Fassung bewenden; bei weiterer Ausbildung wird man natur- 
gemäß auf die erw ähnten mathematischen Grundlagen zurückgreifen müssen. 

Im folgenden werden die Vorgänge in einem dünnen geradlinigen Stab untersucht. Die 
Methode läßt sich auch mit Vorteil auf Bewegungen von Saiten mit variabler Spannung und 
von dehnungslosen Seilen anwenden, worüber ich mehrfach vorgetragen habe und demnächst 
auch im Druck berichten werde. 


2. Problemstellung, Bezeich- 

nungen. In Bild 1 ist schematisch 
ein kurzes Stück eines geraden, Kräfte 
aus homogenem Material bestehen- 

den dünnen Stabes angedeutet, 

dessen Querschnitt F eine stetige 

und genügend oft differenzierbare 

Funktion der jeweils im ruhenden \ 
Raum in Richtung der Stabachse Unstetigkeits-Steik ® 
laufenden Koordinate x ist. Durch Querschnitt F, 
einen bei O— 0 befindlichen Quer- 

schnitt, der sich mit der absoluten 

Geschwindigkeit e in Richtung der 

x-Achse bewegt, wird der Stab in 

zwei Abschnitte 1 und 2 unterteilt. 

In jedem dieser Abschnitte sollen 

die Spannungen o, die Teilchen- 

geschwindigkeiten v, die Dehnun- 

gen e und die spezifische Masse u 

bis an die Stelle 0 — O heran stetig 

und genügend oft differenzierbar 

sein (d. h. so oft, wie es im folgen- OR | | BEE 

den notwendig wird). An der Stelle 

0-0 sollen aber die genannten Bild 1. 

Größen unstetig sein. Sind also 1, 

bzw. I, kleine Längen, die von O—O aus in die Gebiete 1 bzw. 2 gezogen werden, wird 
ferner für alles weitere die Bezeichnungsweise verabredet, daß sämtliche Größen im Ab- 
schnitt 1 den Index 1, im Abschnitt 2 den Index 2 erhalten sollen, und wird ferner 


lim A, = Ay; 
gesetzt, wo A, irgendeine der interessierenden Größen an der durch 1, festgelegten Stelle ist 
(entsprechend für A,), so soll sein: 


In Bild 1 ist dies für o und » schematisch angedeutet. 

An den durch I, bzw. I, gegebenen Stellen seien ferner die Querschnitte a,b, bzw. a, b, 
gelegt, die sich, wie in Bild 1 angedeutet, mit der jeweiligen Teilchengeschwindigkeit be- 
wegen. /, und I, seien von erster Ordnung klein gedacht; sie werden später beliebig klein 
gemacht. Durch die beiden genannten Querschnitte wird ein sich bewegendes Massenelement 
abgegrenzt, das in seinem Innern die Unstetigkeitsstelle enthält. 


3. Die dynamische Grundgleichung. Für dieses Massenelement sei nun zuerst die dyna- 
mische Grundgleichung aufgestellt. Man könnte, wie es meist bei Ableitung der Differential- 
gleichung geschieht, von dem Newtonschen Gesetz ausgehen, würde dann aber an der Un- 
stetigkeitsstelle, an welcher ja ein Spannungssprung besteht, auf Schwierigkeiten z. B. be- 
züglich des Wechselwirkungsgesetzes stoßen, die erst zu beseitigen wären. Viel einfacher 
kommt man zum Ziel, wenn man auf das abgegrenzte Element den Impulssatz anwendet, der 
hiermit in wohl unmittelbar plausibler Weise auch bei Unstetigkeiten als gültig angesehen 
wird. (Die angedeutete Schwierigkeit wird damit allerdings mehr umgangen als beseitigt: 
weiter soll aber hier auf diese‘ Fragen nicht eingegangen werden; die Klärung wäre wohl 
durch einen Grenzübergang herbeizuführen.) 
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Bis auf Größen höherer Ordnung ist der Impuls des abgegrenzten Massenelementes: 


Dies läßt sich ohne weiteres nach der Zeit substantiell differenzieren; für das abgegrenzte 
Massenelement ist also: 


dJ_ 
Hierin sind alle mit !, und !, multiplizierten Größen beschränkt; gehe ich also zur Grenze: 


di 


so bleibt: 
F, Ho’ (+ | 
Offenbar ist: 


so daß schließlich: 


Die auf das abgegrenzte Massenelement in der .-Richtung wirkende Kraft ist bei 
spannungsfreier Oberfläche des Stabes: 


(o,- F) d(o,-F) 


Hierin fallen wieder beim Grenzübergang nach (5) die mit /, bzw. I, multiplizierten Glieder 
fort; aus dem Impulssatz folgt also, da F, sich hier heraushebt: 


Dies ist die erste der zu benutzenden Hauptgleichungen; sie ist von dynamischer Art 


und beruht darauf, daß laut Gl. (7) die zunächst rein geometrisch definierte Unstetigkeitsstelle 
die Eigenschaft einer nicht verschwindenden, wohldefinierten Massenbeschleunigung besitzt. 


4. Die kinematische Hauptgleichung. Eine zweite Hauptgleichung erhalten wir aus der 
Kinematik des Vorganges. 

Die Längenänderung des abgegrenzten Massenelementes erfolgt einerseits nach Gl. (6), 
wobei sofort nach (5) zur Grenze gegangen wird: 


d(,+1,)\ _, 
dt 
Andererseits ist aber die Längenänderung auch durch die Dehnungen gegeben. 


In Übereinstimmung mit der technischen Festigkeitslehre sei die Dehnung & definiert 
durch: 


(8). 


> 
wenn ZL, die ungespannte Länge eines Stabes und Z, die Länge desselben Stabes unter der 
gleichmäßigen Zugspannung o ist. Ist die gespannte Länge gegeben, wie bei !, und I,, so 
ergibt sich hiernach die ungespannte zu 


L 


Sind I,, bzw. I,, die Längen von I, bzw. Il, im ungespannten Zustand, so läßt sich 
schreiben: 


Von höherer Ordnung kleine Glieder sind hierbei bereits fortgelassen. 

Beim Differenzieren von (9) nach t£ ist zu beachten, daß I,, bzw. I,, Konstante sind, die 
nach Wahl von !, bzw. I, festliegen. Geht man dann wieder nach (5) zur Grenze, so fallen 
wieder die mit !, oder !, multiplizierten Glieder fort, und es wird unter Benutzung von (6): 


d(l, _ C— 


_ 
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Mit (8) zusammen ergibt dies: 


C— Un 
1+ Eio 1+ 


als zweite Hauptgleichung. Sie ist rein kinematischer Art. 


(Il) 


3. Das Elastizitätsgesetz. Durch das naturgemäß als bekannt anzusehende jeweilige 
Elastizitätsgesetz 


ist eine dritte Beziehung gegeben, durch welche entweder die Spannungen oder die Dehnungen 
aus (1) oder (II) eliminiert werden können. 


6. Die Gleichungen für die spezifische Masse. Eine weitere Beziehung ergibt sich aus der 
Erhaltung der Masse. Naturgemäß muß sein: 


dt 
also: 


und hierzu tritt schließlich noch eine Gleichung von der Art: 


7. Ein exaktes Ergebnis. Für den vorliegenden Fall kommt man auf folgende Weise zu 
einem ersten übersichtlichen Ergebnis: «, sei die gleichmäßige spezifische Masse des Stabes im 
ungedehnten Zustand. Dann ist, wenn von der Querkontraktion (wie bisher) abgesehen wird: 


Io 


BEER 


= 


da ja 


Ferner sei eine Geschwindigkeit c, dadurch definiert, daß während dt eine Masse 
u, F,-c,-dt von der einen auf die andere Seite der Unstetigkeitsfläche tritt; dies ist aber 
auch F,-(e )dt= F,-(e-v,)dt. Hiermit wird: 


oder mit (11) nach Kürzen mit «,: 


c 


Durch Einsetzen von (12) in (I) uud von (13) in (II) ergibt sich dann: 


und 


Hieraus: 
00 

Durch die Sprunggrößen allein ist also in dem behandelten Fall ohne jede Vernach- 
lässigung die allein für den Massentransport von einer Seite der Sprungstelle zur anderen 
maßgebende Größe c, nach Festlegung einer dazugehörigen spezifischen Masse «, gegeben; 
nicht aber die Fortpflanzungsgeschwindigkeit der Sprungstelle im Raume. 

Gilt das Hookesche Gesetz, so wird mit E als Elastizitätsmodul aus (16): 


E 


Ist der Spannungsdehnungssprung nicht groß und die Kurve o(r) nicht stark gekrümmt, 
so kann oft benutzt werden: 


(iS), 


wobei für den Differentialquotienten rechts ein geeigneter Mittelwert einzusetzen ist. 


8. Die üblichen Vernachlässigungen. In dem hier behandelten Beispiel sind nun folgende 
Vereinfachungen naheliegend und meist zulässig: 
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(19). 


so daß 
c > 
Hiermit wird aus (I) und (II): 


Yu) 


also: 


(23); 
bei kleinen Spannungssprüngen: 
do 
de 
wenn die Kurve o(e) gekrümmt ist. Dabei ist das Vorzeichen von e durch (20) bzw. (21) 
festgelegt. 

Von diesen letzten Gleichungen soll im folgenden hauptsächlich Gebrauch gemacht 
werden. Nach ihnen, d. h. mit den durch (19) ausgedrückten Vernachlässigungen, ist ce, die 
.Fortpflanzungsgeschwindigkeit der Sprungstelle im Raum, eine Materialkonstante, wenn 
zwischen Spannungen und Dehnungen ein linearer Zusammenhang besteht. 


9. Gleichungen für die Geschichte der Sprunggrößen bei Gültigkeit der üblichen Differential- 
gleichung und bei beliebigem Elastizitätsgesetz. Um jetzt Aussagen über das Verhalten der 
Sprunggrößen beim Weiterschreiten der Unstetigkeitsstelle zu erhalten, werde auf die Diffe- 
rentialgleichung zurückgegriffen, die innerhalb der beiden stetigen Gebiete gilt und mit den 
eingeführten Vernachlässigungen (G]. (9)) jeweils lautet: 


(0, F) 0° 
F) 
Hierin ist wie üblich & die Verschiebung; ferner gilt in beiden Indizes: 
_ 98 


dx 


098° 

Das Elastizitätsgesetz wird noch nicht eingeführt. 

Hiermit ergibt sich die Möglichkeit, eine der drei Größen o, &, v aus (21) und (22) zu 
eliminieren, wobei gleichzeitig F als Funktion von x hineingebracht wird. 

Dabei ist folgendes zu beachten: In den Gl. (25) bis (27) erscheinen partielle Differential- 
quotienten der Einzelgrößen; diese existieren mit den betreffenden Differentialquotienten in 
den beiden stetigen Gebieten bis an die Unstetigkeitsstelle heran. 

In den Gl. (21) und (22) kommen aber außer ce und « nur die Sprunggrößen o,, — 0,0; 
— — Vor, und diese existieren nur an der Unstetigkeitsstelle; von ihnen gibt es 
keine partielle Differentialquotienten. Wohl aber haben sie einen Differentialquotienten z. B. 
nach der Zeit, den man erhält, wenn man mit der Unstetigkeitsstelle fortschreitet und die in 
den Sprunggrößen auftretenden Änderungen als Funktion der Zeit auffaßt und differenziert ; 
entsprechende Differentialquotienten gibt es auch für die Einzelgrößen in unmittelbarer Nähe 
der Unstetigkeitsstelle. = 

Es soll nun in Übereinstimmung hiermit Dy diejenige Änderung sein, die eine an der 
Unstetigkeitsstelle oder in ihrer unmittelbaren Nähe existierende Größe y während der Zeit 
Dt erfährt, innerhalb deren die Unstetigkeitsstelle von x, nach ,+Dx,=x2,+e:Dt fort- 


| 
1 0.—0 
und bei Gültigkeit des Hookeschen Gesetzes: 
so daß 
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schreitet. Dann ist hiermit der für e und die Sprunggrößen existierende Differentialquotient 
definiert, also: 


D De 


Derselbe Differentialquotient existiert aber auch für die Einzelgrößen an der Sprung- 
D 
stelle, also z. B. = Ze usw. Öffenbar ist auch z. B.: 
Do, Don 


Dt (0.0 — = 
und umgekehrt; für die anderen Größen entsprechend. 


Für die Einzelgrößen ist aber: 


Dy_(dW 
Dt 


denn hierdurch wird das Mitschreiten mit der Unstetigkeitsstelle zum Ausdruck gebracht; die 
partiellen Differentialquotienten sind jeweils in den stetigen Gebieten in unmittelbarer Nähe 
der Unstetigkeitsfläche zu nehmen. 


Es ist also erlaubt, für die Sprunggrößen an der Unstetigkeitsstelle den Differential- 
quotienten D/Dt zu bilden, in diesem nach (28) die Einzelwerte zu trennen, diese nach (29) 
auszudrücken, sie dann mit Hilfe von (25) bis (27) umzuformen, und schließlich wieder, soweit 
möglich, unter Berücksichtigung von (29) Werte D/Dt der Einzelgrößen zu bilden und zu ent- 
sprechenden Differentialquotienten der Sprunggrößen zusammenfassen. 


Dies wird nun durchgeführt. Wir gehen von Gl. (20) aus und erhalten: 


1D D De 


2 ergibt sich aus (22); die Einsetzung des Wertes in (30) führt mit einfacher Um- 
formung auf: 

D 1 D 


dx 
E 0% 
ersetzen hierin dv/df durch »E und 5 - durch Fr drücken nach Gl. (25) die zweiten Dif- 


d 
ferentialquotienten nach durch aus und beachten, daß F, allein von x, ab- 


Nun bilden wir: 


hängig ist, so daß: 
dinF, DinF, 
DinF,, 


Hiermit wird aus (30) nach Zusammenfassung der beiden Glieder mit Di 


Diese Gleichung läßt sich noch dadurch übersichtlicher machen, daß man auf der rechten 


Seite durch Addition und Subtraktion entweder von oder von 
En dx), 


die Bildung eines weiteren Differentialquotienten D/Dt entweder 
20 10 0 

Von &, —£,, Oder von o,, —0o,, ermöglicht, der dann auf die linke Seite gebracht wird. So 
ergibt sich schließlich das vorläufige Resultat in den zwei gleichwertigen Beziehungen: 
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) 
Dt 
und: 


Auf den rechten Seiten ist eine Zusammenfassung der Differentialquotienten der Einzel- 
größen zu solchen der Sprunggrößen nicht weiter möglich; auf den linken Seiten stehen da- 
gegen nur die Größen 05, — 0,,; &3 —&, und der Querschnitt F, an der Unstetigkeitsstelle. 
Die Zweckmäßigkeit der links gewählten Schreibweise wird sich bei der weiteren Verarbeitung 
der Beziehungen erweisen. 

Man erkennt zunächst: Die Veränderung der Sprunggrößen beim Weiterschreiten der 
Unstetigkeitsstelle ist bei beliebigem Elastizitätsgesetz im allgemeinen nicht allein von ihnen 
selbst und dem Querschnittsverlauf, sondern auch von den Vorgängen in den stetigen Gebieten 
abhängig, aus denen sich die auf der rechten Seite von (34) und (35) stehenden Randwerte 
ergeben. 

Dagegen wird die Geschichte der Sprunggrößen von den Vorgängen in den stetigen 
Gebieten stets dann unabhängig, wenn die rechte Seite der beiden letzten Beziehungen Null 
wird. Auf die Bedingungen hierfür gehen wir nunmehr ein. 


10. Spezialisierung des Elastizitätsgesetzes, insbesondere Hookesches Gesetz; die maß- 


gebenden Gleichungen für die Sprunggrößen als Hauptergebnis.. Wir wählen Gl. (35) als Aus- 
gang und fragen: Wann ist 


Nach Einführung von , worin 0’ die aus dem Elastizitätsgesetz des Stoffes 


zu entnehmende Ableitung von o nach # bedeutet, wird aus (36) mit einfacher Umformung: 


Wir unterscheiden zunächst drei Hauptfälle für das Eintreten von (37). 


a) Bei beliebigem Elastizitätsgesetz und beliebiger Größe des Spannungs- und Dehnungs- 
sprunges sowie der mittleren Spannung und Dehnung stellt Gl. (37) eine Bedingungsgleichung 


zwischen 5) und ($) dar, ist also nur unter ganz speziellen Umständen erfüllt, auf die 
o o 


hier nicht näher eingegangen werden soll. 
b) Ist der Spannungs- und Dehnungssprung bei beliebigem Elastizitätsgesetz „sehr klein“, 
so sind die beiden Klammerwerte in (37) angenähert Null, so daß dann bei nicht zu hohen 


Werten von und 
deutet und wieweit jeweils die Annäherung geht, hängt von der Form der s(e)-Kurve und 
den Besonderheiten des betreffenden Falles ab. 


Für a) und b) wird aus Gl. (35): 
F 
c 


worin €, eine durch Anfangswerte oder dergleichen gegebene Konstante. 
e) Der wichtigste Fall liegt vor, wenn 


Gl. (37) angenähert erfüllt wird. Was „sehr klein“ genau be- 


über den ganzen Spannungs- und Dehnungsbereich, der o,, -o,, und z,, — #,, enthält, wenn 
also in diesem Bereich mit 2 Konstanten a und b: 


also lineare Abhängigkeit herrscht. 
Hiervon wieder liegt der wichtigste Fall vor, wenn a=0, also das Hookesche Gesetz 
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mit E=konst. hinreichend genau erfüllt ist. Dies ist natürlich auch der weitaus wichtigste 
Fall für die gesamte Anwendung der Theorie. 

Unter diesen in der Melırzahl der Fälle zulässigen oder in vielen anderen Fällen infolge 
der sonst auftretenden Schwierigkeiten der Bearbeitung gebotenen Voraussetzung ist Gl. (37) 


immer erfüllt, wie auch die endlichen Werte (5) und (% 


dt 
stetigkeitsstelle ausfallen mögen. 
Da dann auch 


in unmittelbarer Nähe der Un- 


ce — konst., 


so erhält man bei Gültigkeit des Hookeschen Gesetzes mit einer neuen Konstanten (, die 
bemerkenswerte Beziehung: 


Herrscht also etwa zur Zeit i, an der Stelle «,,, beim Querschnitt F,,, der Spannungs- 


sprung 40,, so schreitet dieser mit der Geschwindigkeit e = Vz über den Stab und hat zur 
Zeit f, beim Querschnitt F,, an der Stelle «,,+e(t, —t,) die Größe: 


(43), 


ist also bei geometrischer Ähnlichkeit der Querschnitte irgendeiner gleichgelegenen Länge, 
bei Kreisform also z. B. dem Durchmesser, umgekehrt proportional. Dabei ist der Quer- 
schnittsverlauf, wie aus den Ableitungen hervorgeht, nur der einen Beschränkung unterworfen, 


daß an jeder Stelle = eindeutig ist, die Mantellinie im Längsschnitt also ohne Knick verläuft. 


Aus (20) ergibt sich mit einer neuen Konstanten €, der entsprechende Sprung in der 

Geschwindigkeit zu: 

also die gleiche Abhängigkeit von der Stabforın wie bei o,, — 0o,,. Man beachte jedoch die 
umgekehrte Reihenfolge der Indizes in Gl. (20). 

Eine der Hauptaufgaben dieser Mitteilung erscheint hiermit gelöst. Denn es ist evident, 
daß in solchen Fällen, für die ein unstetiger Vorgang überhaupt als Grenzfall für die Ab- 
schätzung der Verhältnisse brauchbar ist, mit der Kenntnis der Veränderung der unstetigen 
Größen bei beliebigem Querschnittsverlauf schon viel gewonnen ist. Von Sonderfällen ab- 
gesehen, wird man annehmen dürfen, daß z. B. die sonst im Stabe auftretenden Spannungen 
kaum wesentlich größer, sicherlich aber nicht von ganz anderer Größenordnung sein werden, 
als die an der Unstetigkeitsstelle bei ihrem Fortschreiten vorhandenen. 

Einer sinngemäßen Ausdehnung des Resultates auf andere Vorgänge (Torsionsschwin- 
gungen, elektrische und akustische Vorgänge u. a.), die von der gleichen oder einer ähnlichen 
Differentialgleichung erfaßt werden wie das behandelte Problem, steht nichts im Wege. 


11. Ein Beispiel zur Veranschaulichung. An einem einfachen Beispiel sei die Anwendung 
der abgeleiteten Gesetzmäßigkeiten gezeigt. 

Ein Stab der hier behandelten Art bewege sich spannungsfrei mit der für alle Stellen 
gleichen Geschwindigkeit V in Richtung der positiven «-Achse. Zur Zeit t=0 werde der 
Querschnitt am Stabanfang «=0 plötzlich festgehalten. 

Man hat offenbar bei #=0 zur Zeit t=0 einen Geschwindigkeitssprung mit »,, — V und 


v 
dementsprechend nach Gl. (21) einen Dehnungssprung — 


woraus sich mit &,,—0 ergibt: ui und bei Gültigkeit des Hookeschen Gesetzes: 


r 


Diese Spannung entsteht also bei «—0 unmittelbar nach dem plötzlichen Festhalten dieser 
Stelle. Die Störung pflanzt sich mit der Geschwindigkeit ce im Stabe fort, und die Zug- 
spannung im Augenblick des Durchganges der Unstetigkeitsstelle an der Stelle x, mit dem 


(Querschnitt F,, beträgt nach Gl. (42) kann also bei starker Verjüngung 
oi 
des Stabes beträchtliche Werte annehmen. Der Geschwindigkeitssprung ist hierbei nach 
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ze also erhält man unmittelbar hinter der Störungsstelle die Ge- 


Ist der Stab prismatisch, also Vz —1, so wird vı0,=0; der Stabteil links von der 


Gl. (44): 9% = —V: V 
schwindigkeit 


Störungsstelle ist in Ruhe. Bei verjüngtem Stab wird Vr >1, also vıo, negativ; die 


Stabteilchen werden in diesem Fall beim Durchgang der Störung nicht nur zur Ruhe ge- 
bracht, sondern noch nach hinten beschleunigt, entgegengesetzt der Fortschreitungsrichtung 
der Störungsstelle. Bei mit « zunehmendem Querschnitt bleibt vıo, positiv, aber kleiner als V. 
Dies gibt bereits ein ganz anschauliches Bild von dem Vorgang während des ersten Zeit- 
abschnittes, innerhalb dessen die Störungsstelle von «=0 nach dem anderen Stabende läuft. 
Auch das Auftreten der Quadratwurzel in Gl. (43) wird verständlich: Bei verjüngtem Stab 
z. B. geht die auf dem großen Anfangsquerschnitt ausgelöste große Kraftwirkung (Spannung 
mal Fläche) beim Fortschreiten auf kleinere Querschnitte über, die Spannung muß also größer 
werden; gleichzeitig weichen aber wegen des kleinen Querschnittes die verzögerten Massen 
stark nach der negativen «-Richtung aus, so daß die Spannung nicht reziprok dem Querschnitt, 
sondern weniger stark ansteigt. 

Ist das von außen bisher unbeeinflußte Stabende frei, so daß hier keine äußeren Kräfte 
auftreten können, so bleibt der Spannungssprung hier bei seinem Eintreffen erhalten; er 
wandert dann unter Befolgung desselben Gesetzes wieder an die Ausgangsstelle zurück. In- 
zwischen hat sich aber in einem nicht prismatischen Stab unter dem Einfluß der nach dem 
ersten Durchgang der Störung zurückbleibenden Dehnungen und Geschwindigkeiten ein anderer 
Zustand ausgebildet; beim Rücklauf der Störung wissen wir also ohne Integration der Glei- 
chungen für die stetigen Bezirke nichts Bestimmtes über die Einzelgrößen o,,:; 0,0; &s0; Eu; 
%345 do. Die Differenzen bleiben aber entsprechend Gl. (43) und (44) erhalten. 


12. Einige wichtige Einzelpunkte. Folgende Einzelheiten erscheinen noch bemerkenswert: 
f a) Die vereinfachten Gl. (20) und (21) lassen sich anschaulich deuten: Gl. (20) sagt, daß 
durch die Kraft o„,— 0, in der Zeiteinheit eine Masse «-c von der Geschwindigkeit »,, 
auf v,, gebracht wird (Kraft und Masse auf die Flächeneinheit bezogen); Gl. (21) drückt aus, 
daß in der Zeiteinheit eine Länge ce von der Dehnung e,, auf die Dehnung e,, kommt, wobei 
die entsprechende Verlängerung auch durch die Geschwindigkeitsdifferenz der Endflächen 
erhalten wird. 
Aus derartigen einfachen Erwägungen könnte man die Gleichungen wohl auch direkt 
hinschreiben; dabei würde aber die Erkenntnis verlorengehen, daß sie nur angenähert gelten. 


b) Mit Benutzung der Verschiebungsgröße & läßt sich (21) mit (27) leicht umformen zu: 


D 
also 


wobei die Nullpunkte für die Verschiebungen £,, und &, zweckmäßig so gewählt werden, 
daß die Konstante Null ist. Gl. (21) drückt also die Tatsache aus, daß der Zusammenhang 
des Stabes an der Sprungstelle erhalten bleibt, und läßt sich rückwärts hieraus ableiten. 


c) Eine einfache Nachrechnung zeigt, daß die Beschleunigung in F, im allgemeinen eben- 
falls eine Unstetigkeit hat. Doch führt Abwesenheit eines Beschleunigungssprunges nicht 
immer auf im ganzen stetige Verhältnisse. In dem behandelten Beispiel ist z. B. bei prisma- 
tischem Stab die Beschleunigung auf dem ganzen Stab ständig Null, also in F, stetig; dabei 
sind aber die anderen Größen in F, unstetig. 


d) Man kann nach einem Vorschlag von Herrn L. Cremer die Gl. (42) und (44) mit- 
einander multiplizieren und erhält: 


So ergibt sich eine Größe von der Dimension einer Leistung, die konstant ist, wobei F,, in 
der ersten Potenz auftritt. (45) gibt aber im allgemeinen nicht die an der Unstetigkeitsstelle 
bei deren Fortschreiten auf das betreffende Massenelement übertragene Leistung. Diese ist 
vielmehr: 


Lo = > F, Fr . . . . . . . (46). 


2 


| 
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Für diesen Leistungswert gilt, wie unmittelbar plausibel und leicht zu verifizieren, der 
Energiesatz, derart, daß die Arbeit Z,-dt, also die auf das Element von der Länge 1, +1, 
mit 1,0, 1,>0, während dt übertragene Arbeit der an ihm angreifenden Spannkräfte, gleich 
der Erhöhung seiner Energie in derselben Zeit ist, wobei sich diese aus der kinetischen und 
der Formänderungsenergie der beiden Einzelstücke zusammensetzt. Dabei ist im allgemeinen 
L, nicht konstant. 

Dagegen geht (46) dann in (45) über, wenn zwei zusammengehörende Werte von o und v 
Null sind. Es werde z. B. in einem ruhenden Stab der Anfangsquersechnitt plötzlich durch 
eine Druckspannung 9, Mit ®,., in Bewegung gesetzt; jetzt ist beim ersten Wandern der 
Unstetigkeitsstelle über den Stab o,=0 und v,,—=0, also nach (45) und (46): 1, = konst. 

Auf weitere Einzelheiten hierüber soll hier nicht eingegangen werden. 


13. Der mathematische Sachverhalt. Der mathematische Inhalt der vorstehenden Ab- 
scehnitte ist mit den Vernachlässigungen nach (19) der folgende: 


Bei Gültigkeit des Hookeschen Gesetzes in der Form 


o= E d x 
ist die Differentialgleichung des Problems (linear, von 2. Ordnung, hyperbolisch): 
Hierin ist: 
„Il Fe) 
(48), 
also: 
= F(o) - „jonas 
Durch die neuen Variablen 
u=x—ct 
v=xcH4+ct 


sind dann die sogenannten charakteristischen Grundkurven von (47) gegeben, d. h. im £&, «, t- 
Raum die in der x-t-Ebene liegenden Querschnittskurven derjenigen Zylinderflächen, auf 
denen die Charakteristiken liegen; diese selbst ergeben sich erst aus den sonstigen Bedin- 
gungen des jeweiligen Problems, über die hier nicht verfügt wird. 
Wie sich gezeigt hat, existieren jetzt Lösungen von (47) mit folgenden Eigenschaften: 
a) & ist im ganzen Bereich stetig, auch auf den Grundkurven « = konst. und ® — konst., 
so daß in der hier gewählten Schreibweise 


D _ - 


b) Auf den genannten Grundkurven «== konst. existieren Unstetigkeiten in den ersten 
partiellen Differeirtialquotienten, d. h. in den Neigungswinkeln der Tangentialebenen an die 
Integralflächen, so daß also: 


08, 08, 08, 08, 


ce) Als Folge von (51) besteht zwischen den Sprunggrößen folgende Gleichung: 


62), 
ot), dx), \dx 


d) An der Sprungstelle, also auf «—konst. (oder » = konst.) besteht die Beziehung: 


(53). 


e? 


e) Ferner besteht an der Sprungstelle die Beziehung (Energiesatz im mechanischen 
Problem): 
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Stress Systems in an Aeolotropic Rectangular Plate. 
By Hajimu Ökubo in Yonezawa (Japan). 


Verfasser berechnet die Spannungen in einer dünnen äolotropen, recht- 
eckigen Scheibe konstanter Dicke für den ganz allgemeinen Fall, daß die 
vier Kanten durch irgendein System von Zug- und Schubkräften beansprucht 
werden. Er geht so vor, daß er zunächst acht Sonderfälle behandelt, aus denen 
er dann durch Überlagerung die Lösung für den allgemeinen Fall aufbaut. 
Diese enthält im Gegensatz zu früher gegebenen Lösungen nur ein einziges 
System von unendlich vielen Konstanten. Weiter wird, ausgehend von den 
Formeln für die Spannungen in einer unendlich ausgedehnten Scheibe unter 
der Wirkung von Einzelkräften, gezeigt, wie man unter Benutzung der 
gewonnenen Resultate die Spannungen in einer rechteckigen Scheibe be- 
rechnen kann, auf die in einzelnen Punkten Kräfte wirken, deren Rand 
aber spannungsfrei ist. 


$ 1. Introduction. 

The mathematical treatment of the elastic equilibrium of an aeolotropie body in two 
dimensions was carried out by K. Wolf‘) and its further study was developed by the writer?). 
In the preceding paper, the writer introduced the general expression of stress components 
and as an example gave a solution of stresses in an infinitely long strip whose two edges are 
acted upon by some systems of traetions. 

In the present paper, the stresses and displacements in an aeolotropie reetangular plate 
are given for the most general case where the four edges are acted upon by any system of 
traetions. And using the results, the solution is also obtained for the case where isolated 
forces act at the points within the plate and the edges are free from traction. 

From the results of an aeolotropie reetangnlar plate, is readily deduced a satisfactory 
solution of an isotropie reetangular plate under the most general boundars conditions ceited 
above by taking the limiting case when k approaches to unity. 

The problem of determining the stresses and displacements in an isotropie rectangular 
strip was attacked by Filon‘) but the result is not applied without a considerable error 
except the case where the breadth of the strip is small compared with its length. An attempt 
to find the solution for an isolated force acting within an infinitely long strip was made by 
E. Melan*) and a more complete solution was given by R. C. J. Howland?°). The stresses 
in an isotropie reetangular plate whose two edges are acted upon by symmetrically distri- 
buted pressure was caleulated by J. Kuno®). But the solution obtained by him is not con- 
venient for a numerical calceulation since the solution involves four systems of unknown con- 
stants which must be determined by solving four systems of infinite determinants. The 
solution given here involves only one system of unknown constants and is more convenient 
for a numerical caleulation. 


$ 2. Fundamental formulas. 

We consider an aeolotropie rectangular plate of uniform thickness, and take axis of x 
and y in a plane parallel to the faces of the plate and place it so that the boundaries of the 
plate are represented by the lines +a, and y=+b. 

The problem is treated here as one of generalised plane stress, i. e., the stresses and 
displacements considered are means taken with respect to the thickness of the plate. 

The stresses in an aeolotropie body are generally given by’) 

1) K. Wolf: Z. angew. Math. Mech. Bd. 15 (1935), S. 249 bis 254, 

2) H. Ökubo: Seienee Reports of the Töhoku Imperial University (1) 25 (1937), Phil. Mag. April (1939). 

5) L.N. G. Filon: Phil. Trans., A, Bd. 21 (1903), S. 63. 

4) E. Melan: Z. angew. Math. Mech. Bd. 5 (1925), S. 314. 

6) R. ©. J. Howland: Roy. Soc. Proe. A, Bd. 124 (1929), S. 89. 

6 J. Kuno: Phil. Mag. April (1932). 

?) If we denote the stress funetion by Z, the differential equation of X proposed by K. Wolf is 

2 2 
The general solution of the differential equation is 
where f,, f2 are arbitrary functions, and the stress components are given by 
Ay= 
Putting fi’ = Ya, we have the equation (1). 


ch. 
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Y,= 
X,= 


where p, and g, are arbitrary functions which are determined to fit given boundary conditions 
and R and  denote the real part and the imaginary part divided by ö of the complex 
functions respectively and k is the square root of the ratio of Young’s moduli of the x 


and % directions. 
Then the mean displacements « and v are found from 


du 1 


oX,), 
I 
E 


(2), 


where E and o denote Young’s modulus and Poisson’s ratio in the x direction respectively, 


$ 3. Aeolotropic rectangular plate under peripheral tractions. 


We proceed to solve eight special cases from which can be built up the general solution 
satisfying arbitrary boundary conditions of normal stress and shear. Now put 


-+iky)= A, cos" (etiky), 


N 
ylatiy) = B, 


The stress components may be written from (1) 


n x 
X =-k > A„cos B,„«eos — 
a a a 
> Ancos" Y+ Bucos"” Cof 


X,=—k > A,sin Ein” B,sin 


where 
Sinz=—isinixz, 
Case 1. The boundary conditions are 


Z,=0, X,=0, when @=t+ta, 


when y=+b 


The stresses given by (4) will satisfy the boundary conditions 


y= P, cos 


if 


X,=0, when y=+b. 


nay 
a 


a 


(3). 


(4), 


(6), 
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P,„Sin 
EP, &in = b 


But they do not satisfy the boundary ER (5) since 


X,=0, 
X, at the boundary «= + a is an even function of y, we can put 
X =— cos’, 
Then 
P 
where 
b 21.2 2 
\ a a a a b 


Next, putting 


y,(a+iky)= 


j 


3-1 


sin! D;& Sin?” sin 


To satisfy the boundary conditions 


I 


we have 


when z=+a. 
the arbitrary constants C; and D; must be 
Sin, 
Ga 
jra 
za ra 


when z=+a. 


(10), 


(il). 


(12), 


(13). 
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On the edges y= + b, the shear vanishes, but the normal tension is 


j=1 3=1 


which may be expanded in Fourier’s cosine series. 


Rn 
Then 
je 
where 
za 


The stresses (4) superposed with (11) are 
They satisfy the boundary conditions (5), if 
Pn+ Ru=Pn- 
Then P,„ is determined from the following relation. 


P,=P, 


Pu 2 


j=1 


"Enjtin=Pn, 


(14). 


(15), 


(16). 


(m). 


Now P,„ is known, from the equation (9) @; is determined and the stresses (16) are 
known. The displacements u, ® are determined from (2). Then the stresses and displacements 


in any point of the plate under the boundary conditions (5) are all known. 


Case 2. The boundary conditions are 


X,=0, when z=+a, 
X,= Pnsin when y=b, 
Y„=0, Pnsin when y=— b 


We will take the stress system 


NE n kı = 
A, cos > B„cos Coj 


akı 


(18). 


(19). 


_ 
| 
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To satisfy the boundary conditions 


= > pn when y=b, 


pn sin when y=—b 
n—]1 
the arbitrary constants A„, B„ are 
Pn Eoj 
ei) 
— 5 
B,=— 


On the edges = + a, the shear vanishes and the normal tension is 


n—]1l 


which is expanded in Fourier's cosine series. 


where 


Q % 2(— Ir tinpn 
k 
a a >1 
\ b? b? 
By a similar way as Case 1, if we find the stress system X}, Yy, X, which satisfies 
the following boundary conditions 


> X, =0, when +a, 
j=0 . 


y=0, when y=+b 


(23), 


and superpose it to (19), we can find the stress system which satisfies the boundary conditions (18) 
Nextly, we will put 


=D Busin” 


| 

XUM 
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From the equation (1), we have 
k > A, sin Coj > B,„sin 
n=i 
0), 
> A„sin = B„sin Coj 3, 2 
Ser 
> A, cos + > B, cos —, Ein, 
Case 3. The boundary conditions are 
). X,=0, when z=+a, | 
sin =, X,=0, when y=+b | 
The stresses given by (25) will satisfy the boundary conditions 
mu a 
n—| 
if 
nab 
P„Sin 
Aus 
); 
a a a a 
kb (28). 
k Pı Ein 
B 
On the edges x= + a, the normal tension vanishes and the shear is 
kı > 
which may be put 
s 
Pu . . . . . . . . . . . . . . (30), 
n=1 
where 
), 7 kl ah 
2(-1etjjk Sin” Sin — | 
kb ah .nnkb 
Next, putting 
(31), 
= D,;Sın (e+iy) 
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we have 
X =-k > (,;,5in D; Sin” cos 
je 
Yy cos + D, Sin 
j=1 ji= 
b b b 


On the edges x= + a, the normal tension vanishes and the shear is expressed in the 


form X, = \ sin if 
j=1 
9, Sin 
Si Sin — k&üı 
(33). 
Q;k Sin!” 
kb 
IR jrua jra 
On the edges y= + b, the shear vanishes and the normal tension is 
j=1 
which is expressed in the form 
provided with 
=> rj m Q; (35), 
where 
2 
If we determine P, from the following a 
m—| 
then the stress system (25) combined with (32) satisfies the boundary eonditions (26). 
Case 4. The boundary conditions are 
3. X,=0, when 
Y,=0, X, = Pn when y= b, 
n=1 (37). 


Y„=0, = > when y=- b 


n=1 


if 


W 
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The stresses (25) satisfy the following boundary conditions 


X,= Pu €08 when 
(38), 
Y,=®, X,=- > ‚ when y= 
if the arbitrary constants are determined as 
nab 
Pu 
a a a 
(39). 
nn 
B,= 
a a a a 
On the edges x= + a, the normal tension vanishes and the shear is 
X,=k (— 1)" 24 Bu(— 1)” cm 
a a 
“=1 
which is expressed in Fourier’s sine series 
X,=— > 9; sin?” 
je 
where 
2 (— 1)" + 
kb 
k Gof nakb,. 
a a 


By a similar way as employed in Case 8, if we find the stress system satisfying the 
boundary conditions 


2. =0, = Q; sin“, when z=+a, 
x J b ( 4) 
=b, when y= +b 
and superpose it to the stress system (25) eombined with (39), we can obtain the 'stress 
system which satisfies the boundary conditions (37). 
We will find another type of stress system by putting 
y,(a+iky)=i > A, cos 


(42). 


=i > cos 


The stresses are 
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nnky .nnz,.. nn 
A „sin —— Sin —— — \ B„sin 
a ne a a 


Y,= > En + > 
n a n a a 


\' naky n nay 
a Coj + \ B, 


Coj 


A,„cos 
— a 


Case 5. The boundary conditions are 
X,=0, 


Xy 


Pn sin a X; y> 


nn 
er 
Y, Pn sın a ’ 


The stresses (43) satisfy the boundary conditions 


N Pu sin X,=0, 
— a 
s=]- 

u 


Pu sin when 


kb 
_ 


B„= 
kb 7 
Sin nn 


When z= + a, the normal tension vanishes but the shear is 


X,=k 1" B„(— 1)" &o 


which is 


=] 


k Sin 


similar way as employed in the next case satisfying 


If we find a stress system b a 
the boundary conditions 


| | 
| 
Xy=0, when y=—h | | 
| 
Yy= when y= 
n—=] | 
if 
z 
A,„= | 
a ‚nakb b 
a a 
where ® 
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= > cos when @= +a, 


when y= +b 


and superpose it to (43), we have a stress system which satisfies the boundary conditions (44). 
Case 6. The boundary conditions are 


%,=0, whenz=+ a, 


Y,„=0, X,= Pn €0S ‚ when y=+b 


n—1 
The stresses given by (43) satisfy the following boundary conditions 


= 
> P, eos when y= + b 


Sin 


a 


— P, Sin nnkb 
a 
B,= 
nakb nnb „. nnkb nzb 


Sin Sin 


a a a a 


The normal tension vanishes on the edges x — + a, but the shear is 


or 


n—l 


which is expressed in the form 


Jay 
(;c0s 


jel 


Neun, 


n—] 
where 


‚.nnb 
2(- tit' Sin Sin 


7 b 


a 


Next, putting 


| 
if 
Sin” 
then 
n?k? 
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we have 


sin?’ + D; Sin 7, 


To satisfy the boundary conditions 


\' 
; cos 


X, when +a 


the eonstants C,; and D; must be 


Sin k Sin kb Coj 


jra jaa 
a 
Coj kb Sin Sind; 


On the edges y= + b, the normal tension vanishes and shear is 


Koh 


which is 


and 


where 


a 
Sin? 


The stress system (43) combined with (55) satisfies the boundary conditions (49), if 


Pa 


Lastly, we will put 


Ansin” 


tiymi > B,sin 


j=1 
j=] j=1l 
:& ; 
X, > + > D;&oj cos, 
j=1 j=] 
. . . (7). 
| Ä Ä | 
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The corresponding stresses are 


B„cos 
( 


nnc_.. nıaky 
X, Sin F 
3 = a a 
\ nze_. nnky nnx 
A, cos Sin \ B„cos 
a a 


Case 7. The boundary conditions are 


X,=0, X,=0, when z=+a, 
— 


\ Pu €0S ‚X,=0, wheny=-—b 
a 


The stress system (62) satisfies the boundary conditions 


Pn €0S X,=®, when y= b, 


X,=0, when y=-b 


Au 


nnkb_ nnb nab' 
Sin Sin 
a a a a 


- kp. Coj 
B,= 


kC&of Sin 


nakb, nab „ nakb. .nab 
Ein —- —& Coj 
a a a 


a 


On the edges «= + a, the shear vanishes and the normal tension is 


| 


n—] n—| 
which is 
\ Q;sın 
where 
nnkb.. nnb nıkb 
- m — Sin — 


a 


b kb 5 ; 
- So Sin 
a a a ( 
k? 
a? b? 
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The stresses which satisfy the following boundary conditions 


. 
Q; sın when «= d, 


Yy=0, X=0, when y=+b 


(67), 


are found by a similar way as ÜUase 3. If we superpose them to the stresses given by (62), 
we can find a stress system which satisfies the boundary conditions (63). 


Case 8. The boundary conditions are 


,=%; when «= + 


. 
Y=0, X= when y= + 


The stresses (62) satisfy the boundary conditions 


y NNnE 
> Pusin when y= + 


‚nnab 
Pu Ein 
nakb_ nab nah‘ 
Sin —— 


a a a a 


nakb 
nab _. nnkb nab 


— Sin - Ein Eoj 


a a a a 


Pu Sin 
By 


On the edges £—= + a, the shear vanishes and the normal tension is 


A„(- Sin v 1)" Sin 


Jay 
Q; sin 


a 


’ 


which is 


where 


0; = 
ab? Coj - Sin "Ein —— &of 


a 


nakb, nab 
9m 
a 


By the same way as Case 7, if we find a stress system satisfying the boundary con- 
ditions (67), and superposing it to the stress system (62) combined with (70), we can find the 
stress system which satisfies the boundary conditions (68). 


$ 4. Solution under general peripheral tractions. 


Firstly, we will assume that the upper and lower edges of the plate are acted upon by 
any distributions of normal force and shear, while the lateral edges are free from tractions. 
The boundary conditions of the case are 


= 
| 
if 
Au= 
J 
| 
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Yıına 
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; 
Y„= cos + B, sin 
a a 


n=u 


X,= C„ os — + D,„ sin 
a a 


X,= N cos + D’, sin 
u a a 


n=] 


when y=b, 


when y=—-b 


which is readily transformed into 


na“ 
B,) sin 22 — B,) sin a 


when y=b, 


n 
(On + cos 


\ (A„+ 4/,) cos — A,) cos 


a 


u 


when y=-—b 


D;,) sin 


Hence the boundary condition (72) is the combination of eight boundary conditions (5), 
(18), (26), (37), (44), (49), (63) and (68). The stress system which is built up by adding the 
solutions of the eight cases is the solution satisfying the boundary condition (72). By a 
similar way, if we find a solution of the case where the upper and lower edges of the plate 
are free from tractions and the lateral edges are acted upon by any distribution of normal 
force and shear, and adding both solutions together, we can find the solution under general 
peripheral tractions. In the boundary condition (72), tlıe terms of €, and €, must be equal 
since the peripheral traetions of the plate are in equilibrium. But when the plate is under 
the peripheral tractions acting on the four edges, the tractions on opposite edges need 
not be in equilibrium respectively, if only the resultant is in equilibrium. 


$ 5. Forces operate within the plate. 
We will firstly find the stresses in a plate unbounded in all direetions corresponding 
to an isolated force P which acts at the origin in the positive direction of x. 
Now, putting 


| 
| 
(RB), 
| 
. | J 
1 
n=]1 
N 
n—1 . (3). 
1 
1 
n—1i 
1e 
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| (74) 


we have 


X,=— 


Ax Bx 
. 
y 

The resultant of the tractions is found by adding the stresses on the cavity bounded 
by a eircle with its centre at the origin, since the resultant is independent of the shape of 
the cavity. The equation (75) shows that the component of the resultant in the direetion 
of y vamishes, and the moment of the traetions about the centre of cavity also vanishes, and 
the component in the direetion of x is 27(k A+ B). 


Hence we have 

To make the displacement v one-valued 

A+(l+o)kB=0. 

Pkil+o) 
27 (k? —1)’ 
Similarly, when an isolated force Q acts in the positive direction of y, the corresponding 
stresses are 


Yy= 


then, we have 
A= 


B= 


x’ +k?yp ’ 
B' x 


X,=- 
Y,= 


X,= 


where 
—1)’ 
When several forces act at the points within the plate and are in equilibrium, and 
the force JP;, Q;) acts at the point (w;, %;), the corresponding stresses in the plate are 


4’ 


By 


Aa) 


where 


i 2a(k®—1) ’ 

If we find the stress system, by the way explained in $4, in a rectangular plate under 
the peripheral traetions, which are equal and opposite with the stresses (80) on tlıe edges 
of the plate, and superposing it to the stress system (80), we can find the stresses in the 
plate which is acted upon by forces within the plate but is free from tractions on all edges. 

In eonelusion the writer wishes to express his hearty thanks to Prof. Iwao Kobayashi 
and to Prof. Seiichi Higucehi of the Töhoku Imperial University for their kind guidance 


in this study. 216 
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Über unvollständige Anger-W ebersche Funktionen'). 
Von Peter und Elfriede Brauer in München. 
(Aus den Physikalischen Instituten der Techn. Hochschule München und der Bergakademie Freiberg i. Sa.) 


Die unvollständigen Angerschen und Weberschen Funktionen werden 
dureh ihre Integraldarstellung eingeführt. Eine Reihendarstellung wird für 
sie gegeben. Für einen wichtigen Sonderfallwerden Tafeln im IntervallO <o< x 
berechnet, die sich aber — wie abschließend gezeigt wird — mit Hilfe ein- 
facher Formeln auch bei beliebigem reellen Argument verwenden lassen. 


1. Einleitung. Versucht man die durch ihre natürliche Gleichung ?) °) 


bestimmte Kurve in rechtwinkligen Koordinaten « und »® darzustellen, so ergibt sich folgende 
Parameterdarstellung 


Setzt man in den Integralen für die obere Grenze o—=n ein, so bemerkt man, daß man — 


bis auf den fehlenden Faktor S — im Falle (2a) die Angersche Funktion J,(r) und im 


Falle (2b) bis auf den Faktor — — die Webersche Funktion E,(r) vor sich hat®). Bekannt- 


lich gehen ihrerseits die Ange vo Funktionen für ganzzahliges p in Besselsche Funk- 
tionen über. — Wir wollen in (2) setzen u=u/,(o) und v=v/,(o). Es gilt also 


und 


2. Reihenentwicklung. Wir wollen für die Funktionen (2) Reihenentwicklungen angeben. 
Wir erläutern unser Vorgehen an der Funktion «. Es ist 


u=| cos (rsino— po)do=\ cos (r sin o)cospodo-+ \sin (r sin o)sinpodo, 


Hierin lassen sich cos (r sin o) und sin (rsino) nach Jacobi in Fouriersche Reihen ent- 
wickeln: 
eos (r sino)=J,(r)+2 an (r)cos2no 


sin (r sin 0) — 2 (r)sin@n+1)o. 


Da diese im Intervall O<o=<27 gleichmäßig konvergieren, tun es auch diejenigen Reihen, 
die aus ihnen durch gliedweise Multiplikation mit dem konstanten, d.h. von » unabhängigen 
Faktor cospo bzw. sin po hervorgehen. Wir können also gliedweise integrieren: 


o 


u=J,(r)\eospodo+2 (r) |cos2n ocospodo 
0 


+2 
0 


‘) Bezeichnung nach Vorschlag von Herrn Prof. Fritz Emde. 
2) Es bedeutet K die Krümmung, © die Bogenlänge der Kurve, p und r Konstante. 


3) Diese Gleiehung trat bei Verfolgung einer technischen Aufgabe auf. Siehe hierzu P. Brauer, Straße 
Bd. 7 (1940), S. 352 bis 357. Die Behandlung der Aufgabe, bessere Übergangsbögen im Straßenbau auf Grund kine- 
matischer Überlegungen anzustreben, wurde von Wa. Ostwald (Straße Bd.7 (1940), S. 144 bis 145) angeregt. 


4) Siehe z.B. R. Weyrich: Zylinderfunktionen, Leipzig 1937, S. 87 oder ©. N. Watson: A Treatise on the 
Theory of Bessel Funetions, Cambridge 1922, S. 308. 
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Die Integrale können nach Teilen integriert werden, so daß wir erhalten 


J,(r)sinpo +2 > "J,n(r) (sin 2nocospo — 37608 2no sin po) 
1 


0 
In analoger Weise ergibt sich 


1 


v 
Setzen wir hierin o=x, so erhalten wir 


x 


» ha J,(r) —2  +Sinpr . (a), 


Diese Reihen müssen nach (3) mit den Weberschen Reihen für 7 J, bzw. - a 
p + (2? — — (4? — (2? p°) (6? — 


+sinpa 


r 


übereinstimmen. Zum Beweise setzen wir in (5) für J die zugehörige Potenzreihe 


> 
m)! 


ein und ordnen die Summenglieder in (5) innerhalb einer Klammer dann nach Potenzen 
von r. Dabei nn sich die Glieder mit r“ (unter 0 des in (5) vor der Klammer 


bzw. — zu 
) p 


1 


wenn u=2m+2n gerade ist. Ist u—2m+(2n+1) ungerade, so haben wir (unter 
Weglassung des Faktors sinpr bzw. (1 + cos px)) 


stehenden Faktors — 


“—] 


1 
2 (-1)2 


r" 2 - — 
1 


5) H. F. Weber: Ann. Phys. (III) Bd. 8 (1879), S. 407 bis S. 444 siehe S. 419 und S. 420. 
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Diese Ausdrücke stellen aber die Partialbruchzerlegungen der entsprechenden Glieder in (6) 


(2? pP?) (4? pP?) (u gerade) 


p*) 


4 (u ungerade) 


dar, wie man durch eine kleine Rechnung zeigen kann. 


In ähnlicher Weise kann man versuchen, die Reihen (4) für Bruchteile von x zu be- 
rechnen, doch erhält man keine so einfachen zur numerischen Rechnung geeigneten Ausdrücke. 


3. Tafel. Setzt man besonders p=r, so erhält man den für die Anwendungen wichtigen 
Fall, bei dem die Anfangskrümmung der Kurve, deren rechtwinklige Koordinaten « und v 


sind, verschwindet. p hat dann die Bedeutung des mit — 4 multiplizierten Winkels zwischen 


den Kurventangenten in den Punkten o=0 und o=n°). Für diesen wichtigen Fall wird 
die Tafel gegeben für eine Reihe von Werten von p zwischen 0 und 0,5. Ihre Berechnung 
geschah mittels tabellarischer Integration. Das Intervall O<o=<7 wurde dazu in 20 gleiche 


Teile der Breite == geteilt und für jedes o=n 55 y, der Integrand in (2), mit Hilfe einer 


siebenstelligen sin- und cos-Tafel’) berechnet. Das Differenzenschema wurde bis zu den 
6ten Differenzen ausgedehnt, wozu die Integrandenfunktion auf beiden Seiten des Intervalles 
um je drei Intervalle erweitert werden mußte. Zur Berechnung der Teilintegrale mit den 


Grenzen ?n und n+ 95 sowie derjenigen mit den Grenzen (2 und (2n-+3) 
wurde die Formel°) 


1 1 


20 
verwendet, zur Berechnung der „Anschlußintegrale“ \ bzw, | hingegen die Formel’) 
19 


20 


un 


_ 
5 » 


70 2 60402 


4+4, 191 4+ 
2 


Nach diesem Verfahren ergaben sich die Integralwerte mit 8 zählenden Stellen nach dem 
Komma. Unsere Tafel wurde auf 6 Stellen abgerundet. Eine Kontrolle ließ sich bequem 
durchführen, da die Berechnung der Endwerte (s=r) auch mit den Reihen (6) durchgeführt 
werden konnte'®). Die festgestellten Abweichungen betrugen im Höchstfalle weniger als 
3 Einheiten der 8. Stelle nach dem Komma. Deshalb wurde bei den Teilintegralen, bei 
deren Abrundung eine um den jeweils möglichen Fehler geänderte 8. Stelle die 6. beeinflußt 
hätte, diese unter- oder überstrichen. Eine unterstrichene 6. Stelle ist also schlimmstenfalls 
um 1 Einheit zu groß, eine überstrichene wäre zu klein. 


4. Erweiterung für beliebiges Argument. Es sei darauf hingewiesen, daß Funktionswerte 
jeder der beiden Funktionen für beliebiges reelles Argument sich linear durch die Funktions- 
werte beider Funktionen im ersten Duanten (d.h. O<o=<7) ausdrücken lassen. 

R 6) Siehe P. Brauer a.a.0. S. 354 Gl. (20a). 

?) Opus Palatinum. Sinus- und Cosinus-Tafeln. Herausgegeben von Dr. W. Jordan, Hannover und Leipzig 1913. 

s) Siehe z.B. ©. Runge und H. König: Vorlesungen über numerisches Rechnen, Berlin 1924, S. 239 Gl. (1). 

») ©. Runge und H. König: a.a.0. S. 241 Gl. (2). 


10) Hierbei wurden die Wstelligen Tafeln, herausgegeben vom Reichsamt für Landesaufnahme unter wissen- 
schaftlicher Leitung von Prof. Dr. J. Peters, Berlin 1919 und 1922, verwendet. 
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Bild 1. 
Bild 1. Zur Bestimmung von Funktionswerten für beliebiges reelles 
Argument. Die Ziffern 1 bis 3 bezeichnen den jeweiligen Nullpunkt des 
ein- bis dreimal gespiegelten Koordinatensystems, dessen v®-Achse die 
Verbindungslinie von Z mit dem Nullpunkt ist. 


Bild 2. Unvollständige Angersche Funktion. 


Bild 3. Unvollständige’Webersche Funktion. Darstellung der Funktions- 
werte mit umgekehrtem Vorzeichen. 


Bild 2. 


ir 


Man gelangt dazu durch folgende elementare Überlegung. Die Kurve © (u, ») = 0 (Bild 1) 
setzt sich in den Intervallen 2nz <o<(2n+1)r aus gleichen Kurvenstücken, in den Inter- 
vallen aus hierzu spiegelbildlich gleichen Kurvenstücken zusammen, 
wie man sofort aus der Periodizität der natürlichen Gl. (1) schließen kann. Man gelangt, 
wie man am einfachsten aus Bild 1 abliest, durch Spiegelung des Koordinatensystems an 
der Kurvennormalen nz (Bild 1) zu einer neuen relativen Lage von Kurvenstück und Koordinaten- 
system, welche die zu den Argumentwerten des 2ten Duanten (r<So=<27) gehörigen Funktions- 
werte « und » als Koordinaten der Kurvenpunkte des ursprünglichen Kurvenstückes ablesen 
läßt. Für den nächstfolgenden Duanten 2r<o=<37) gelangt man durch Spiegelung des 
ersten Spiegelbildes an der Kurvennormalen », im Anfangspunkt des Kurvenstückes zur 
richtigen Lage des Koordinatensystems. Durch sukzessive Befolgung dieser Vorschrift er- 
hält man für die Funktionswerte « und » im „ten Duanten (mn —N)r<o=<nz) auf elemen- 
tarem analytischen Wege folgende einfache Ausdrücke 


(n gerade) 

u(o)= —u(o’)cos(n + (o)— sinn +V)pz, (n ungerade) 
bzw. 

r(0)= u(o)sinnpr + (n gerade) 

v()= tv — v*}cos(n+l)pr — v*, (n ungerade). 
Hierin ist 0’ der ein ganzzahliges Vielfaches von = überschreitende Betrag der Bogenlänge, 
also nam +0'=no. wu(o’) und v(o’) sind demnach die der Tafel zu entnehmenden Funktions- 
werte im iten Duanten. »* ist die Konstante 


Natürlich kann man auch, wie Bild 1 lelırt, sehr schnell durch Konstruktion mit Zirkel und 
Lineal zu den weiteren Funktionswerten gelangen, welches Verfahren zur Herstellung der 
Bilder 2 und 3 benutzt wurde. 

Unser Dank gebührt den Herren Prof. Fritz Emde (Stuttgart) und H. Poeverlein 
und Dr. H. Voit (München) für nützliche Ratschläge und den Herren Prof. R. Tomaschek 
(München) und Prof. G. Grüß (Freiberg i. Sa.) für die freundliche Erlaubnis zur Benutzung 
von Rechenmaschinen. 257 
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Über das Photo eines Photos. 


Von Ulrich Graf in Danzig. 

Herrn Prof. Dr. E. Salkowski zum 60. Geburtstag gewidmet. 

Wird ein Lichtbild eines räumlichen Motives von einer beliebigen Stelle 
aus unter willkürlicher Aufnahmerichtung photographiert, so ist das ent- 
standene Bild im allgemeinen nicht mehr identisch mit einer Original- 
aufnahme des räumlichen Gegenstandes; dies letzte ist vielmehr nur für 
zwei ausgezeichnete Aufnahmerichtungen an jeder Stelle des Raumes der 
Fall. Das erste dieser beiden Photos gibt die Reproduktion des Lichtbildes, 
das zweite eine neue Originalaufnahme. Diese beiden erstrebten Bilder 
können durch Umzerrung aus einer willkürlichen Aufnahme gewonnen 
werden. Aus jedem Photo eines Photos läßt sich daher für jeden Aufnahme- 
ort außer der Reproduktion des Lichtbildes noch eine einzige wohlbestimmte 
„gefälschte Originalaufnahme“ des räumlichen Motives herstellen. 


I. Die Photographie eines Lichtbildes. 

Bei einem Lichtbildervortrag wurde aus dem Zuschauerraum das gezeigte Lichtbild photo- 
graphiert. Es entstand die Frage, ob das so gewonnene „Photo eines Photos“ noch identisch 
sein könne mit einer Originalaufnahme des im Lichtbild vorgeführten räumlichen Gegen- 
standes bzw. — da die perspektiven Transformationen ja nicht die Gruppeneigenschaft be- 
sitzen — welche besonderen Aufnahmebedingungen eingehalten werden müßten, um ein solches 
Photo eines Photos von einer Originalaufnahme ununterscheidbar zu machen. Diese Frage, 
deren Beantwortung die Beurteilung von photographischen Fälschungen u. ä. gestattet, ist im 
folgenden untersucht. Gegeben sei also ein Lichtbild eines räumlichen Gebildes. Aus einem 
Photo dieses Lichtbildes soll eine neue Originalaufnahme des räumlichen Motives gewonnen 
werden. 

Eine Photographie möge dabei durch eine exakte zentralperspektive Abbildung ersetzt 
und weiterhin zunächst keine Rücksicht darauf genommen werden, ob das zentralperspektive 
Bild in solcher Nähe des Bildhauptpunktes entsteht, daß seine Realisierung mit einem der 
üblichen Photogeräte möglich ist. Auf diese letzte Frage wird zum Schluß eingegangen. 

Das erste Photo [P] des räumlichen Gegenstandes [®] entstehe durch Zentralprojektion 
vom Zentrum O0 auf die Bildebene = (Bildweite d, Hauptpunkt H) (Bild 1). Jedes weitere 
Photo dieses Bildes kann, wenn überhaupt, nur identisch sein mit einer Originalaufnahme 
des räumlichen Gegenstandes, die von O als Zentrum aus gemacht ist, denn eine Aufnahme 
des räumlichen Motives von einem anderen Zentrum aus würde Einzelheiten der räumlichen 
Gestaltung aufweisen, die auf dem Bilde x und damit auch auf jedem Photo von ihm nicht 
zu sehen sind. Es bleibt somit zu untersuchen, welche Aufnahmen des z-Bildes identisch 
sind mit einer Aufnahme dieses Bildes von O selbst aus. Diese letzten Aufnahmen liefern 


ia, da sie als ebener Schnitt des von O ausgehenden Strahlenbündels zu dem räumlichen 


Gegenstand entstehen, stets Bilder, die als Originalaufnahmen des Raumgegenstandes an- 


gesehen werden dürfen. 
Bild 2 (rechts). 


Bild 1 (unten). 
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II. Die ausgezeichneten Aufnahmerichtungen. 


Alle zentralperspektiven Bilder, die mit einer Aufnahme von O aus identisch sind, liefert 
das bekannte Parallelogrammgesetz der Perspektive (Bild 2). Ist die erste Auf- 
nahme durch das Zentrum O, den Hauptfluchtpunkt F, den Hauptspurpunkt S und den Haupt- 
verschwindungspunkt V gekennzeichnet, so denkt man sich diese vier Eckpunkte des Parallelo- 
gramms OFSV als Gelenkpunkte. Jede kinematisch zulässige Lage 0'’F’SV, die durch Be- 
wegung des Gelenkparallelogramms über SV als fester Grundstrecke hervorgeht, führt auf 
eine neue Aufnahme (Zentrum 0’, Hauptfluchtpunkt F’, Hauptspurpunkt S, Hauptverschwin- 
dungspunkt V), die mit der Ausgangsaufnahme (O0 FS V) übereinstimmt. Zur analytischen Er- 
fassung dieses Sachverhaltes wird ein «yz-Achsenkreuz benutzt, wie es Bild 1 zeigt: die 
yz-Ebene ist die Bildebene x, das Zentrum O hat die Koordinaten =d, y=z=0. Jede 
perspektive Aufnahme wird durch ihren „Feldvektor“ 5 gekennzeichnet; der Aufpunkt (x, y, 2) 
ist das Zentrum der Perspektive, die Richtung von ‘% gibt den Hauptstrahl und die Länge 
von % die Bildweite an. % ist somit der Stellungsvektor der Verschwindungsebene. Bild 2 
zeigt die Lageverhältnisse in einer Ebene durch die x-Achse mit den Koordinaten & und 
w=yy’+z’. (Die Vektoren $ und — % in demselben Aufpunkt liefern perspektive Bilder, 
die nur die bekannte Vertauschung von links und rechts, oben und unten aufweisen und da- 
her als identisch angesehen werden dürfen.) 

Alle Aufnahmen im Hauptzentrum O werden durch die Vektoren %, gekennzeichnet, die 
stets die feste Bildweitenlänge |%,|=f haben mögen. Sind «a,, £,, y. die Winkel eines solchen 
Vektors %, mit den drei Koordinatenachsen, so ist demnach 


cos’c,+ cos’, + cos’y,—1 


Die mit den Aufnahmen { 5,} übereinstimmenden Aufnahmen {75} an irgendeiner anderen 
Stelle x, 9, liefern insgesamt ein Vektorfeld, das rotationssymmetrisch bezüglich der «-Achse 


ist. Die Feldlinie für jeden Ausgangsvektor %, ist der in der Ebene 


gelegene Halbkreis durch 0, dessen Mittelpunkt in x liegt und für den %, Tangentenvektor 
ist. Wird der Drehwinkel « als Parameter gewählt, so lassen sich die Feldlinien in der Form 


(08 a, + sin 


=, , (C08 a, + sin u 
sin’ a, ) 


darstellen. Der Feldvektor im Punkt (x, y, 2) ergibt sich als Tangentenvektor der Feldlinie 
unter Beachtung der elementargeometrisch an Bild 2 leicht zu bestätigenden Tatsache, daß 
31:19] 


ist, bei geeigneter Orientierung zu 
f 1 


A. Im Falle cos’ cos?’ = cos? y,—=0 7) wird die in (3) gegebene 


Darstellung für y und z unbestimmt; in jedem Punkt (x,»,2) hat dann nach (5) der Auf- 
nahmevektor die Gestalt 


d 
008%, 
sin a, 
XUM 


mn 


wr Zu 
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(Bei entgegengesetzter Orientierung würde man die gleichwertigen Aufnahmevektoren 
x -i erhalten.) Dieser triviale Sonderfall liefert also alle Aufnahmen, bei denen die 


Hauptstrahlen parallel zum Hauptstralil des vorgelegten Photos sind; die neu erzeugten Photos 
sind einfache Reproduktionen des Ausgangsphotos. Das Vektorfeld (#’) ist ein Potentialfeld, 
denn es ist 


mi 


Bild 3a zeigt das Vektorfeld in einer Ebene durch die «-Achse. Der zu jedem Vektor % ge- 
hörige Grundvektor %, ist = — fi. 


9=Wonst Bild 3a (links). 


Bild 3b (rechts). 


d [A296 3e] 


B. Ist «+0; 7, so ergibt die Elimination von a,, ß,, y. aus (1), (3) und (5) den Feldvektor 


f 


Die Winkel a, ß, y dieser Aufnahmerichtung gegen die drei Achsen sind also durch 


cos a: cos B:cosy=(ar— 
) 
cos’ a+ cos’ p+cos’y=1 | 
gegeben. Ferner gilt 


sina:sina,=x:d, cos ß:cos y= ß,:C0S 


Die Vektoren % bilden ein komplex-lamellares Feld, denn es ist 


3=4:grad 
mit 
2x 
und 
Ver+y+z+ d’? — 4d? x? 


Der Maßstabsfaktor A bewirkt, daß der Betrag von % zur Bildweite einer Perspektive wird, 
die mit der zugehörigen Aufnahme {75,} von O aus identisch ist. 
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Die Feldlinien sind‘ die in (3) dargestellten Kreise, die Orthogonalflächen ® = const 
sind, der geometrischen Deutung entsprechend, die Kugeln des elliptischen Kugelbüschels mit 
den Grundpunkten (Nullkugeln) = +d, y=z=0. Bild 3b zeigt das Vektorfeld in einer 
Ebene durch die x-Achse. Der zu jedem Feldvektor % gehörende Grundvektor %, ist der 
Tangentenvektor der Feldlinie im Punkte 0. Er ergibt sich zu 


Die durch ihn festgelegte Aufnahme des räumlichen Motives stimmt mit der Aufnahme {5} 
des vorgelegten Photos überein. Die Aufnahme {7%} erfolgt dabei unter den Winkeln «,, 
Po; y. gegen die Achsen des Bezugssystems, gegeben durch 


a,:008 B,:008 =(@ +? +2? — d’):2dy:2dz, | (13) 
cos’ a, + cos? B, + y, = 


Der dureh %, bestimmte Feldlinien-(Halb)kreis hat den Mittelpunkt V mit den Koordinaten 


cosa,.coß, _y 


(14) 
und den Halbmesser 


An jeder Stelle des Raumes gibt es . 
zwei Aufnahmerichtungen, unter denen 
sich ein Photo oder Projektionsbild eines 
räumlichen Motives derart photographieren 
läßt, daf diese neuen Aufnahmen mit Ori- 
ginalaufnahmen des räumlichen Motives 
übereinstimmen. Unter Benutzung des in 
Bild 1 eingeführten Koordinatensystems 
sind die beiden Aufnahmen im Punkte 
(x, y,2) durch die Vektoren und [(6) 
und (9] gekennzeichnet (Bild %. Die 
ersten Aufnahmelagen ‘5 ergeben als tri- 
vialen Fall die einfachen Reproduktionen 
(ähnliche Bilder) des vorgeführlen Licht- 
bildes; in ihrer Gesamtheit bilden die 
Vektoren '%' ein Potentialfeld (Bild 3a). 
Die zweiten Aufnahmelagen % ergeben vom 
Ausgangsbild abweichende Photos, die auch 
unmittelbare Aufnahmen des räumlichen Bild 4 
Motives, und zwar von O aus unter den 
Winkeln a,, ßs, Y. (13) gegen die drei Achsen, sein könnten. Die Vektoren 5 bilden in 
ihrer Gesamtheit ein komplex-lamellares Feld (Bild 3b). Die Bildweiten || und |%| sind 
so gewählt, daß die zugehörigen Originalaufnahmen von 0 aus die Bildweite f besitzen. 


III. Die Umzerrung einer willkürlihen Aufnahme. 


Eine andere Wahl der Bildweiten, durch die sich stets Potentialfelder erreichen ließen, 
würde nur eine ähnliche Vergrößerung oder Verkleinerung des Photos bedeuten; solche „ge- 
fälschten Originalaufnahmen“ lassen sich daher mit Kammern jeder Bildweite machen, und 


2 
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die erforderliche Einstellschärfe läßt sich stets erzielen. Bei Aufnahmen mit starker Neigung 
gegen die Ebene des Projektionsbildes würde die Scheimpfugbedingung zu beachten sein. 
Jedoch scheitern im allgemeinen solche Aufnahmen zunächst praktisch an dem Umstand, daß 
dann das Bild weit entfernt vom Hauptpunkt liegt und auf der Platte unerreichbar wird. 


Die Hauptpunkte der vier Aufnahmen seien bzw. Zi, Z, Z,Z 
(Bild 4); der Hauptpunkt H des Projektionsbildes hat seine vier Bilder in Z/, P,, P’', P. Diese 
vier Bildpunkte liegen jeweils auf denjenigen Geraden, die auf ihren Plattenebenen durch 
die Ebene (2) (Zeichenebene in Bild 4) ausgeschnitten werden. Eine kurze Rechnung ergibt 
nach (10) und (13) 


2 


(16), 


wobei x der Winkel des zum Aufpunkt 0’ führenden Radiusvektors gegen die &- Achse ist. 
Somit ist bei den beiden Aufnahmepaaren {%,}, {%, und {750}; {75} der Streckenunterschied / 
zwischen den Hauptpunkten Z und den Bildern P des Punktes H, d. h. des Hauptpunktes 
des vorgelegten Projektionsbildes, beidemal der gleiche, und zwar ist er proportional dem 
Abstand des Aufnahmeortes (x, y,2) vom Hauptstrahl „=z=0. Wenn daher das entstehende 
Bild auf dem begrenzten Plattenrechteck überhaupt Platz finden soll, so lassen sich die vor- 
geschriebenen Aufnahmen {75} und {7%} nur in der unmittelbaren Umgebung dieses Haupt- 
strahles ausführen. Das aber würde eine starke Einschränkung der Möglichkeit bedeuten, 
mit den üblichen Geräten ein Photo so zu photographieren, daß das neue Bild ein Original- 
photo des räumlichen Motives sein kann. 


Um das Plattenrechteck voll auszunutzen, wird man das Photogerät dem Projektions- 
lichtbild zuwenden, d. h. den neuen Aufnahmevektor ® in 0° durch den Hauptpunkt H dieses 
Lichtbildes gehen lassen (Bild 5). Dann liegt ® auf dem Radiusvektor nach 0’, hat also bei 
der gewählten Orientierung die Form 


17), 


wobei » die Bildweite des benutzten Gerätes bedeutet. Die so erzeugte Aufnahme {3} läßt sich 
nun wiederum in die beiden geforderten Aufnahmen {%} und {%} umzerren, die ihrerseits 
mit den Originalaufnahmen {%,} und {75,3 übereinstimmen. (In Bild5 sind die Vektoren %, 
und % in der Gegenrichtung derer in Bild 4 eingetragen, was jedoch nach dem Vorstehenden 


nur die Vertauschung von links und rechts, oben und unten bedeutet.) 
Aus den drei Plattenebenen der Aufnahmen {2}, {5°} und {7} schneidet die Ebene (2) 


drei Geraden aus, die das Dreieck BD B’ bestimmen. Für die Dreieckswinkel »’ = X DB'’B 
und = X DBB’ ergibt sich elementargeometrisch nach Bild 5 


und nach (10) und (17) 
+2+d) 


2 2 
tg — vY go= 


Aus der Aufnahme {3} lassen sich daher durch Umprojektion auf geneigte Bildebenen die 
beiden geforderten Aufnahmen {7%} und {7%} herstellen, am bequemsten in einem der photo- 
grammetrischen Entzerrungsgeräte. Bild 6 zeigt noch einmal die beiden Einstell-Lagen, wobei 


t 
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der Bauart der Geräte gemäß die Punkte 4’, B’, 0", 0’ 
bzw. A,B,C, 0’ als Gelenkpunkte zu denken sind, die 
eine kinematische Verschiebung dieser Parallelogramme 
gestatten und so einen Freiheitsgrad für die Scharf- 
einstellung gewährleisten. In Erfüllung der Scheimpflug- 
bedingung muß die Hauptebene des Objektives 0’ die 


Bild >. 


Gerade 0’ B’ bzw. 0'B enthalten. Die erste Lage Bild 6a liefert dann die Entzerrung zu einer 
einfachen Reproduktion des Ausgangsphotos, wie es z. B. bei dem „Wiederaufrichten“ lotrechter 
Kanten an Architekturaufnahmen aus der Ballon- oder Froschperspektive gebräuchlich ist. 
Die zweite Lage Bild 6b liefert ein Photo, das mit einern Originalphoto des räumlichen 
Motives, aufgenommen von O aus unter den Winkeln a,, f,, Y. (13) gegen die drei Achsen, 
übereinstimmt. Zwischen den Winkeln «a,, a, ®, & besteht unabhängig von x, 9, z, d außer 
(18) nach (10) und (13) die Bindung 


tra 


2 


(20). 


sin 20’)=sina: + sin?’ ctg 


Der geometrische Ort aller Punkte (x, y, 2), von denen -aus die Aufnahmen {7%} unter 
dem gleichen Winkel a zu machen wären, ist nach (10) die Fläche : 


2°). 
Sie entsteht durch Rotation der gleichseitigen Hyperbel 


um die x-Achse; alle diese erzeugenden Hyperbeln bilden ein Büschel, das von der Hyperbel 
und dem Geradenpaar <= aufgespannt wird. 

Der geometrische Ort aller Punkte (x, y, 2), von denen aus Aufnahmen für den gleichen 
Winkel a, gemacht werden können, ist nach (13) der Torus 


(y’-+ 2°). 


Alle diese Flächen entstehen durch Rotation der Feldlinien (3) um die x-Achse. 


Wird die Aufnahme {®} ohne Verkantung ausgeführt, also derart, daß zwei Kanten des 
Plattenrechtecks waagerecht liegen, so ist für die Entzerrung die Kenntnis des Winkels v er- 
forderlich, den die Schnittgerade der Ebene (2) und der Plattenebene mit den Höhenlinien 
der letzten bildet (Bild 7a). (In Bild 2 bis 6 ist stets die Ebene (2) als Zeichenebene ge- 
wählt.) Für » erhält man (Bild 7b) 


(21). 
Wird ein Photo oder Projektionsbild von irgendeiner Stelle (x, y, 2) aus mit zu- 
gewandter Photokammer photographiert, also derart, daß der Hauptstrahl des Gerätes 
durch den Hauptpunkt H des Lichtbildes geht, so läßt sich diese Aufnahme {8} in die ge- 
forderten Aufnahmen (75, und {75} umzerren. Die beiden Einstell-Lagen und die Neigungs- 
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winkel w' und & der Ebenen sind durch (19) und Bild 6 festgelegt. Die erste Lage ergibt 
demnach die einfache Reproduktion des Projektionsbildes, die zweite Lage ein Photo, das 
mit einer Originalaufnahme des räumlichen Motives von O aus unter den Winkeln a,, Ps; 
y. (13) gegen die Achsen übereinstimmt. 


Bild 7a. 


Bild Sa (unten). 


Bild 8b (Mitte). 
Bild Se (rechts). 


Die Bilder 8a, b, c zeigen als Beispiel solcher Photos eines Photos die Krypta in der Kirche 
Worms-Hochheim. (Aus reproduktionstechnischen Gründen sind die Bilder als Durchzeichnungen wieder- 
gegeben.) Das vorgelegte Photo der Krypta hatte die Bildweite d=8,5 cm, rekonstruiert nach be- 
kannten elementaren Gesetzen der Zentralperspektive. Die Aufnahme {®} dieses Photos für 2=z 
=212em, y=0, v—=W°, |B|=30 em ist in Bild 8a wiedergegeben. Die 
Umzerrung des Bildes 8a nach den obenstehenden Regeln ergab einerseits die Aufnahme {%°), also die 
Reproduktion des Ausgangsphotos (Bild 8b) mit |%°|= 21,2 cm, andererseits die Aufnahme {#} (Bild 8e) 
mit 5 =21,2 em, =50°%, a y=5". Dieses Bild 8e stimmt mit einer Originalaufnahme 
der Krypta von O aus mit = 24°, =W°', [=d=8,5 em überein. 


Aus jedem Photo eines Photos läßt sich somit außer der Reproduktion des Ausgangs- 
bildes noch eine einzige wohlbestimmte Originalaufnahme des räumlichen Motives gewinnen. 
Ihr wahrer Aufnahmeort ist stets derselbe wie der des vorgelegten Lichtbildes, die Aufnahme- 
richtung jedoch hängt von der Wahl der Stelle ab, von der aus das Lichtbild photographiert 
wird. Werden von einem räumlichen Motiv, etwa einem Bauwerk, zwei verschiedene ge- 
eignete Lichtbilder vorgeführt und (von einer Stelle oder auch von verschiedenen aus) photo- 
graphiert, so lassen sich demnach diese beiden „Photos eines Photos“ so umbilden, daß sie 
eine bequeme stereoskopische Auswertung in einem der photogrammetrischen Geräte ge- 
statten. 256 
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KLEINE MITTEILUNGEN 


Bemerkungen zu einigen bestimmten 
Integralen. In dem bekannten Tafelwerk „Nou- 
velles Tables d’Integrales Definies“ von Bierens 
de Haan, Leyden 1867, sind einige bestimmte 
Integrale, die man in der Elastizitätstheorie der 
unendlich ausgedehnten Halbscheibe benötigt, 
falsch angegeben. Es handelt sich um folgende 
Riechtigstellungen: Es muß heißen: 


0 (Hauptwert) für p 


\ dr 
)@—-)yı- a 
0 


2 


statt yp für |p| >1 in table 12 (1), 
! gür p<i 
\ In’ 2 
/ r? 


statt -aln2p für p/>1 in table 119 (1). 
1 


Für \ m dzistintable 119 (2) nur der 


0 


p 1 angegeben. Man braucht aber gerade den 


Wert für >1, der sich zu 2) 


ergibt. Diese und ähnlich gebaute Integrale lassen 
sich sofort angeben, wenn man p als komplexe 
Variable ansieht und beachtet, daß: 


| 
\/(@)-In R(\/(e) In (1 — dx) 
0 
gilt, sowie: 

Ip 


| 


Beispiel: (1 x?) 


pr+1 
00 


p 

1 

2aresine + laresin (7,2) 
pl 


— aresin dp 


Inp+ + | ı) 


am 


Für reelles p erhält man bei p <1 den oben 
angegebenen Wert, während bei p >1 


2 2 > 


folgt. Derselbe Wert ergibt sich auch, wenn man 
im Reellen rechnet und den Cauchyschen Haupt- 
wert bildet. 


München. Gerhard Schubert. 266 


Zur Beschleunigungskonstruktion an 
der Heusinger-Steuerung. Die Geschwindig- 
keits- und Beschleunigungskonstruktion an der 
Heusinger-Steuerung wird im Sehrifttum?) mit 
ähnlichen Punktreihen ?) und fiktiven Konstrukti- 
onen ?) durchgeführt. 

Für die Geschwindigkeiten hat der Verfasser schon 
ein unmittelbares Verfahren ohne ähnliche Punkt- 
reihen veröffentlicht *°), die Beschleunigungskon- 
struktion aber wird in der vorliegenden Arbeit 
behandelt, wobei der Verfasser hier, wie auch in 
anderen Fällen ®°), den Plan relativer Normal- 
beschleunigungen benutzt ®). 


1) M.Grübler: Getriebelehre. Berlin 1917, Seite 81 bis 87. 


2) R. Beyer: Technische Kinematik. Leipzig 1931, Seite 
281 bis 283. 


3) F. Wittenbauer: Graphische Dynamik, Berlin 1923, 
Seite 517 bis 523. 


N. Rosenauer: Eine unmittelbare Gescehwindigkeits- 
konstruktion der Heusinger-Stenerung für Lokomotiven. 
Acta Universitatis Latviensis. Mech. IN 15, 1936. 


5) N. Rosenauer: Über die Gesehwindigkeitskonstruk- 
tion kinematischer Ketten. Z. angew. Math. Mech. Bd. 17 
(1937), S. 173 bis 176. 


6, N. Rosenauer: Die Beschleunigungskonstruktionen 
kinematischer Ketten mit Hilfe von Plänen relativer 
Normalbeschleunigungen. Acta Universitatis Latviensis. 
Mech. II N 4, 1937. 


?) N. Rosenauer: Beschleunigungskonstruktionen kine- 
matischer Ketten mit Hilfe von Plänen relativer Normal- 
beschleunigungen. Masch.-Bau, Reuleaux - Mitteilungen 
Bd. 6 (1938), S. 543 bis 546. 


N. Rosenauer: Ein Beitrag zur Beschleunigungs- 
konstruktion acht- und mehrgliedriger Gelenkketten. Acta 
Universitatis Latviensis. Mech. II N 14 (1939). 


®) N. Rosenauer: Zur Beschleunigungskonstruktion 
acht- und mehrgliedriger Gelenkketten. Masch.-Bau, Ge- 
triebetechnik Bd. 7 (1939) S. 557 bis 560. 
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Bild 3 und 4. 


digkeiten V,; und V,, des Punktes Ä 
und des Kulissenzentrums P,, kon- 
struiert, darauf wird die gedrehte 


Geschwindigkeit des Hilfspunktes 
Q (Schnittpunkt der Geraden AP,, mit 
dem Gliede 5) gefunden, der als Punkt 
des Gliedes 6 aufgefaßt wird. Weiter 
kann aber die Konstruktion in dieser 
Richtung nicht fortgesetzt werden, da 
der Drehpol des 6. Gliedes nicht bekannt 
ist. Daher wird wieder von der Ge- 
schwindigkeit des Radpunktes H aus- 
gegangen, die gedrehte Geschwindigkeit 


7 des Punktes G konstruiert und auf 
den Punkt F übertragen. Darauf wird 


die gedrehte Geschwindigkeit V, des 
Hilfspunktes R (Schnittpunkt des 9. 
Gliedes mit der Normalen zur Gleit- 
richtung des Schiebers im Punkte E) 
gefunden, der als Punkt des Gliedes 7 
aufgefaßt wird. 


Die Heusinger-Steuerung besteht aus 
12 Gliedern, die nach Bild 1 bezeichnet 
sind. 

Die Abmessungen und Lagen ein- 
zelner Glieder sind so angenommen, 
daß die graphische Konstruktion deut- 
licher hervortritt. 

Es wird angenommen, das Rad 
(Glied 1) drehe sich gleichförmig mit 
bekannter Geschwindigkeit gegen den 
Rahmen, (Glied 

Zu ermitteln sind die Gleitbeschleu- 
nigung des Schiebers und die Be- 
schleunigungen der Gelenkpunkte. 

Um die Beschleunigungen zu kon- 
struieren, muß man vorher die Ge- 
schwindigkeiten kennen. Dieselben 
sind in Bild 1 ermittelt worden, wobei 
einfachheitshalber die gedrehten Ge- 
schwindigkeiten der Radpunkte gleich 
dem Abstande vom Radzentrum ange- 


nommen sind V;=JOundV„=H0O, 
dann sind auch die Beschleunigungen 
und a„=HO. 

Die Geschwindigkeitskonstruktion 
soll hier ganz kurz beschrieben werden 
(genauer s. Fußnoten 4 und 5). 

Von der bekannten Geschwindigkeit 
des Punktes / ausgehend, werden 
üblicherweise die gedrehten Geschwin- 


Bild 1 und 2. 
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Die erhaltenen Geschwindigkeiten Vg und V, ge- 
ben die Möglichkeit zum Gelenk C, das die Glieder 6 
und 7 verbindet, überzugehen und dessen gedrehte 


Geschwindigkeit V. üblicherweise zu konstruieren. 
Nachdem diese Geschwindigkeit erhalten ist, findet 
man auch alle >; gedrehten Geschwindigkeiten 


des Gliedes 6 gegen XI. 
In Bild 1 sind auch nachträglich die Drehpole 

und der Glieder«6 und 7 als Schnitt- 


Letztere in der Bewegung 


punkte der Richtung von v. mit dem Gliede 5 und 
der Normalen n —n gefunden. 

Zur Ermittlung der fernse_ nötigen Coriolis- Be- 


schleunigung findet man „konstruiert die re- 


o 
lative Geschwindigkeit V„ nach der Formel 
74XU 73XU 
Va Va, 


vorbiistnt den Punkt M mit dem Eigene von 


zieht durch den Endpunkt vonV„ eine 
parallele Gerade, die die Gerade P,, A im Punkte L 
trifft. Der Abschnitt ‚AL ist laut “Proportion 


Mi A 
die halbe Coriolis-Beschleunigung, da 


2 


Ferner bildet man die Geschwindigkeitsdifferenzen 
und findet mit der Halbkreiskonstruktion die rela- 
tiven Normalbeschleunigungen JK,.HG,FD,.BT, 
der Glieder JA,HG,FD,BT uud KM, A, B. "BC, 
DC’, C,E der Geraden KM, AB, BC, „pe, CE. 


Ebenso konstruiert man mit Normal- 
beschleunigung AN des Punktes A in der Relativ- 
bewegung 4 gegen 3 

Der Plan relativer Normalbeschleunigungen ist 
deutlichkeitshalber in einem zweimal größeren 
Maßstabe (Bild 2) folgendermaßen gezeichnet: 
ijo—=2JO.ik=2JK,, km=2KM,, Akma»AKMA, 
ab=2AB, cd=2CiD, 
df=2D,F, /h=2G,H, da die Punkte / und y im 
Plane zusammenfallen und endlich 5t=2BT.,. 

Die Grundeigenschaft dieses Planes besteht darin, 
daß die Endpunkte der Beschleunigungs- 
vektoren zweier Punkte eines Gliedes auf 
einer Geradenliegen, die zur Verbindungs- 
geraden dieser Punkte senkrecht ist !",!t), 

Die Konstruktion der Beschleunigungen am er- 
haltenen Plane wird folgendermaßen durchgeführt: 
die Geraden durch den Punkt o_Lik und durch 
mk gezogen ergeben den Endpunkt von 
Eine Gerade durch den letzteren L ka und_eine 
durch m L ma ergeben den Endpunkt von jaXI, 

Von demselben trägt man laut Formel 


die Coriolis-Beschleunigung j,=4 AL, die Normal- 
beschleunigung der Relativbewegung j,=2AN 
ab und zieht 1 AP,, die Gerade a—ca, auf der 
sich der Endpunkt von j{X!! befindet. 


10), N. Rosenauer: Die Beschleunigungskonstruktionen 
Kur u Ketten mit Hilfe von Plänen relativer 
Normalbeschleunigungen. Acta Universitatis Latviensis. 
Mech. II N 4, 1937. 

11) N. Rosenauer:Beschleunigungskonstruktionen kine- 
matischer Ketten mit Hilfe von Plänen relativer Normal- 
et. Masch.-Bau, Reuleaux - Mitteilungen, 
Bd. 6 (1938), S.543 bis 546. 


Durch den Punkt ? zieht man _LBT die Gerade 
auf der sich der Endpunkt _von befindet. 

Vom Punkte A trägt man j„—=ho,—=2HOabund 
zieht durch o, eine Gerade 1_HG bis zum Schnitt- 
punkte mit der Geraden, die durch den Punkt /(g) 
parallel zur Gleitrichtung des Gliedes X gezogen ist. 
Der Abschnitt der letzteren ergibt direkt j7= jg, durch 
deren Endpunkt die Gerade ö&— ö_ LFD gezogen 
wird, auf der sich der Endpunkt von ja befindet. 

Die letzte dieser Geraden e— €, auf der sich der 
der Endpunkt von j. befindet, wird direkt durch 
den Punkt e parallel zur Gleitrichtung des Schiebers 
(Glied 8) gezogen. 

Weiter wird nach dem Verfahren des _. 
sers!?) zwischen den Geraden a—a und 
beliebig die Gerade 4’B’_L AB gezogen, auf der- 

BE 
AB AB 
gefunden und mit dem Schnittpunkte U’ von a—a 
und 3—2 verbunden. 

Ebenso wird zwischen den Geraden ö—ö und 
e — & beliebig die Gerade D’E’_LDE gezogen, auf 
derselben der Punkt C’” laut der Proportion 

DE' DE 
gefunden und mit dem Schnittpunkte U” von d— ö 
und e—e verbunden. 

Der Schnittpunkt C, der Geraden C’U’ und C”U” 
ergibt den Endpunkt der Beschleunigung des Ge- 
lenkes C, so daß = je ist. 

Durch den Punkt C, zieht man noch A,|C’ 4’, 
D,E,\D’E’ und erhält die Beschleunigungen der 
übrigen Gelenkpunkte: ad, =), 

dD,=ja und eE,=|).. 

Es muß noch bemerkt werden, daß die Be- 
schleunigungen in Bild 2in einem zweimal größeren 
Maßstabe erhalten sind. 

Anmerkung. Bei gewissen Abmessungen einzel- 
ner Glieder und in gewissen Lagen kann es vor- 
kommen, daß die Geschwindigkeiten der Punkte 
der Glieder DE oder AC nahezu parallel sind. 
Es sind dann die Geschwindigkeitsdifferenzen sehr 
klein und die relativen Normalbeschleunigungen 
verschwinden. Das bedeutet, daß einige Punkte 
im Plane relativer Normalbeschleunigungen zu- 
sammenfallen. Dieser Umstand bietet jedoch kein 
Hindernis zur Beschleunigungskonstruktion nach 
diesem Verfahren. 


Die Bilder 3 und 4 sen einen Fall, 


selben der Punkt C’ laut Proportion 


wo die 


Geschwindigkeiten Va, und nahezu pa- 
rallel sind. 

Es fallen deshalb die Punkte c,d,e,/,g_im Plane 
zusammen, die Beschleunigungen aber je, ja, je und 
jr weisen einen beträchtlichen Unterschied auf 
und der Beschleunigungspol I’ des Gliedes DE 
ergibt sich auf der Normalen n—n. 

Es sind zur Kontrolle in Bild 3 die Beschleunigun- 
gen je, jd. je und js in entsprechendem Maßstabe 
aufgetragen und man kann sich davon überzeugen, 
daß die Endpunkte von jg,je und j. auf einer 
Geraden liegen, und daß die Projektionen dieser 
Beschleunigungen auf die Gerade ED gleich sind. 


Ebenso sind auch die Projektionen von jq und 
jr auf FD einander gleich. 
Riga. 


12), N. Rosenauer: Zur Beschleunigungskonstruktion 
acht- und Masech.-Bau, Ge- 
triebetechnik Bd. 7 (1939), S. 557 bis 560 


N. Rosenauer. 238. 
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